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Presently  available  laser  velocimeter  (LV)  electronic  Bignal  pro- 
cessing techniques  are  often  inadequate  for  detection  of  light 
scattered  by  small  scattering  particles  which  are  required  for  fol- 
lowing the  fluid  motions.  Example  situationr  include  low  pressure 
and/or  high  velocity  measurements.  In  other  situations  detection 
of  larger  scattering  particles  is  difficult  due  to  limited  system 
sensitivity.  Photon  counting  techniques  offer  improved  system 
sensitivity  by  allowing  velocity  measurements  to  be  made  even  when  J 
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^there  are  Insufficient  signal  photons  available  to  define  the  clas- 
sical scattering  signal.  Such  techniques  are  thus  applicable  when 
the  presently  used  classical  techniques  are  not. 

The  objective  of  this  contract  was  the  theoretical  and  experimental 
development  of  a new  type  of  photon-counting  processor  for  low- 
level  dual-scatter  LV  signal  detection.  The  contract  was  conducted 
in  three  separately  funded  phases:  the  first  was  a mathematical 

feasibility  and  preliminary  hardware  design  study  for  a photon- 
counting processor;  the  second  was  a more  detailed  study  of  the 
statistical  nature  of  low-level  LV  signals  and  their  effects  on 
low-level  signal  processing:  and  the  third  phaBe  was  the  detailed 
design,  construction,  and  demonstration  of  a new  photon  counting 
signal  processor  for  photon  resolved  laser  velocimeter  signals. 

The  effort  under  this  contract  has  resulted  in  several  significant 
conclusions  both  theoretical  ana  experimental.  , 
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PREFACE 

This  , technical  report  covers  work  performed  under 
Contract  No.  F49600-74-C-0016  for  the  Arnold  Engineering  « 
Development  Center,  Arnold  Air  Force  Station,  Tennessee. 

The  work  was  performed  between  May  3,  1974  and  October  31, 
1975  by  the  Atlanta,  Georgia  Office  of  Science  Applications, 
Inc.  The  program  was  under  the  direction  of  Captain  M.  L. 
Guiou,  Arnold  Engineering  Development  Center.  ' The  Program 
Element  Number  was  65807F. 

The  reproducibles  used  in  the  reproduction  of  this 
report  were  supplied  by  the  author,  W.  T.  Mayo,  Jr.,  Science 
Applications,  Inc. 
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SUMMARY 

Presently  available  laser  velocimeter  (LV)  electronic 
signal  processing  techniques  are  often  inadequate  for  detec- 
tion of  light  scattered  by  small  scattering  particles  which 
are  required  for  following  the  fluid  motions.  Example  situ- 
ations include  low  pressure  and/or  high  velocity  measurements. 
In  other  situations  detection  of  larger  scattering  particles 
is  difficult  due  to  limited  system  sensitivity.  Photon  count- 
ing techniques  offer  improved  system  sensitivity  by  allowing 
velocity  measurements  to  be  made  even  when  there  are  insuffi- 
cient signal  photons  available  to  define  the  classical  scat- 
tering signal.  Such  techniques  are  thus  applicable  when  the 
presently  used  classical  techniques  are  not. 

The  objective  of  this  contract  was  the  theoretical  and 
experimental  development  of  a new  type  of  photon-counting 
processor  for  low-level  dual-scatter  LV  signal  detection. 

The  contract  was  conducted  in  three  separately  funded  phases: 
the  first  was  a mathematical  feasibility  and  preliminary  hard- 
ware design  study  for  a photon-counting  processor;  the  second 
was  a mere  detailed  study  of  the  statistical  nature  of  low- 
level  LV  signals  and  their  effects  on  low-level  signal 
processing;  and  the  third  phase  was  the  detailed  design,  con- 
struction, and  demonstration  of  a new  photon  counting  signal 
processor  for  photon  resolved  laser  velocimeter  signals.  The 
effort  under  this  contract  has  resulted  in  several  significant 
conclusions  both  theoretical  and  experimental.  These  are 
summarized  below. 

Theoretical  Poisson  models  have  been  formulated  for  dual- 
acatter  LV  signals  with  sufficient  generality  to  include  all 
signal  ’'egimes  from  classical  signals  with  Gaussian  noise 
down  to  photon-resolved  signals  which  require  photon  counting 
techniques . 
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Analysis  of  conditional  signal  models  has  provided 
instantaneous  signal-to-noi se  ratio  and  other  statistical 
parameter  equations  not  previously  available  in  the  LDV 
literature. 

Analysis  of  Unconditional  models  has  indicated  the  prac- 
tical feasibility  of  improving  detection  sensitivity  over 
burst-counter  techniques  by  a factor  between  100  and  1000  by 
using  photon  counting  techniques. 

A detailed  analysis,  including  explicit  mention  of 
approximations  required,  of  the  expected  value  of  direct  pho- 
ton cori jlation  of  LV  signals  from  turbulent  flow  is  provided. 
This  both  supports  and  points  out  limitations  of  previous 
results  given  by  Pike  and  others. 

A new  Dual-Correlate-and-Subtract  technique  for  new  fre- 
quency measurement  which  utilizes  only  two  delay  values  and 
one  real-time  high  speed  multiplier  without  clipping  has  been 
theoretically  shown  to  be  feasible.  This  feasibility  is 
based  on  evaluation  of  the  expectation  for  a system  model  to 
show  convergence  to  mean  frequency  estimate,  and  on  rn  rms 
variability  analysis  which  indicates  practical  measurement 
times  on  the  order  of  0.5  seconds  in  transonic  flow. 

An  electronically  switchable,  multipurpose,  1-rgh-speed 
electronic  signal  processor  with  the  potential  sp  jd  for  10 
sec  time  resolution  and  10  nsec  real  time  4 bit  multiplication 
has  been  designed  and  constructed  using  slow  ECL  integrated 
circuits  on  special-purpose  100  MHz  wirewrap  panel.  The  sys- 
tem implements  the  dual-correlate-and-subtract  technique. 

It  also  has  the  potential  of  standard  photon  correlation 
without  clipping  by  computer  controlled  scanning  of  the 
delays.  The  computer  interface  has  not  yet  been  constructed, 
however.  If  the  computer  interface  is  developed,  the  system 
will  also  be  capable  of  measuring  multipoint  statistics  of 
any  high-speed  events  (up  to  100  MHz) . 


4 


AE  rtC-T  R'76'81 


Experimental  measurements  have  been  performed  which 
verify  the  operation  of  the  photon  processor  with  100  MHz 
input  events  and  40  MHz  system  clock  operating  speed.  A 
sequential  photon  correlation  of  a low-level  sinusoidally 
modulated  optical  source  has  been  made  with  25  nsec  time 
resolution.  The  mean  velocity  of  a low-speed  jet  was  meas- 
ured with  the  dual  correlate  of  subtract  mode  under  condi- 
tions for  which  no  burst  counter  data  could  be  obtained. 

Photon  counting  experimental  work  requires  much  more 
care  in  selection  and  use  of  photomultiplier  tube  components 
than  does  burst-counter  work.  However,  single  photoelectron 
pulse  height  statistics  and  dynode  collection  efficiency  can 
easily  reduce  the  effective  signal  power  to  noise  power  ratio 
by  a factor  of  4 for  two  tubes  with  the  same  quantum  efficiency. 
This  can  significantly  effect  a burst  counter  system  sensitiv- 
ity . 

The  statistical  distribution  of  the  amplitudes  and  rates 
of  occurrence  of  classical  bursts  has  been  shown  to  be  cen- 
tral in  the  problem  of  specifying  or  predicting  the  data 
rates  and  errors  from  any  type  of  LV  signal  processor.  Dif- 
ferential and  cumulative  rate/amplitude  distributions  have 
been  formulated  and  analyzed  theoretically  and  have  been  meas- 
ured experimentally  for  an  argon  backseat ter  LV  system.  The 
results  indicate  that,  for  the  data  obtained,  the  smaller 
aerosols  contribute  more  to  the  photon  correlation  accumula- 
tor than  the  larger  ones.  For  the  data  measured,  there  would 
have  been  available  less  than  300  signals  per  second  adequate 
in  magnitude  to  produce  burst  counter  data  from  scatterers 
larger  than  0.7  micron  in  diameter  while  there  would  have  been 
over  100,000  signals  per  second  producing  photon  resolved 
signals  from  0.2  - 0.3  micron  diameter  particles. 

We  have  demonstrated  both  theoretically  and  experimen- 
tally that  there  are  LV  signals  that  require  the  use  of  photon 
counting  techniques;  and  we  have  delivered  a new  type  of 
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photon  counting  instrument  which  will  potentially  operate 
to  10  nsec  time  resolution  with  full  4 bit  multiplication. 

Further  effort  with  computer  controlled  data  collec- 
tion is  required  to  explore  the  potential  of  the  new  system 
and  fully  define  its  limitations. 
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1 . 0 INTRODUCTION 


1.1  THE  PROBLEM 


Presently  available  laser  velocimeter  (uV)  electronic 
signal  processing  techniques  are  often  inadequate  for  detec- 
tion of  light  scattered  by  small  scattering  particles  which 
are  required  for  following  the  fluid  motions.  Example 
situations  include  low  pressure  and/or  high  velocity  meas- 
urements. In  other  situations  detection  of  larger  scattering 
particles  is  difficult  due  to  limited  system  sensitivity. 
Photon  counting  techniques  offer  improved  Bystem  sensitivity 
by  allowing  velocity  measurements  to  be  made  even  when  there 
are  insufficient  signal  photons  available  to  define  the 
classical  scattering  signal.  Such  techniques  are  thuB  applic- 
able when  the  presently  used  classical  techniques  are  not. 

The  objective  of  this  contract  was  the  theoretical  and 
experimental  development  of  a new  type  of  photon-counting 
processor  for  low-level  dual-scatter  LV  signal  detection. 

The  contract  was  conducted  in  three  separately  funded 
phases:  the  first  was  a mathematical  feasibility  and 

preliminary  hardware  design  study  for  a photon-counting 
processor;  the  second  was  a more  detailed  study  of  the  sta- 
tistical nature  of  low-level  LV  signals  and  their  effects 
on  low-level  signal  processing;  and  the  third  phase  was  the 
detailed  design,  construction,  and  demonstration  of  a new 
photon  counting  signal  processor  for  photon  resolved  laser 
velocimeter  signals. 

During  the  conduct  of  the  contract  a complementary 
theoretical  effort  with  the  NASA  Langley  Research  Center 
was  in  progress  (NAS1-13140  and  follow-on  NASAl-13737 ) . 

The  original  feasibility  mathematics  developed  in  the  first 
phase  of  this  contract  have  been  refined  in  the  NASA  effort. 
Portions  of  this  report  incorporate  some  of  the  NASA  results 
for  conciseness. 
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1 . 2 BACKGROUND 
1.2.1  Signal  Modeling 

Earlier  modeling  efforts  have  treated  LV  signals  for 
which  the  noise  could  be  considered  as  additive  independent, 
stationary,  and  Gaussian  [1,2,3].  This  is  the  limiting  case 
of  stationary  Poisson  shot  noise  which  occurs  for  visible 
light  photodetection  when  a steady  light  source  such  as  a 
heterodyne  reference  beam  [1]  or  high  background  light 
level  [2]  dominates  the  signal.  In  a recent  simulation  of 
low-level  dual  scatter  signals,  the  accuracy  of  the  noise 
model  was  extended  by  treating  it  as  a non  stationary  GausBian 
process  whose  variance  is  proportional  to  the  incident 
optical  power.* 

When  we  discuss  "noise"  in  LV  signal  detection  we  are 
usually  referring  to  the  variation  of  the  electronically 
detected  signal  with  respect  to  a scaled  version  of  the 
classical  optical  (power)  signal  incident  on  the  PMT.  In  a 
general  analysis  the  classical  LV  signals  are  also  random 
processes  due  to  the  random  amplitudes  and  the  arrival  times 
of  the  signal  bursts.  Mayo  [1]  has  treated  these  signals  sb 
a Poisson  process  for  steady  flow,  and  Durrani  [3]  and  George 
[5]  have  treated  them  for  the  turbulent  flow  case  as  Gaussian 
in  the  limit  of  high  particle  number  density.  A new  book  by 
Snyder  [6]  treats  generalised  Poisson  processes  in  great 
detail.  In  particular  he  treats  "doubly  stochastic  Poisson 
processes."  These  are  inhomogeneous  ( n on  stationary)  Poisson 
processes  for  which  the  rate  function  is  a random  process. 
Such  a description  is  appropriate  either  for  the  classical 
LV  signal  bursts  with  the  random  turbulent  flow  affecting 
the  rate  of  burst  occurrence  or  for  the  single  photo-electron 
pulses  from  the  PMT  with  the  random  classical  bursts  as  the 


*Modified  version  of  noise  model  described  in  [4]. 
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rate  function.  Clearly,  when  taken  from  the  turbulence  to 
the  photo-electron  pulses,  a dual-scatter  LV  signal  is  a 
"triply  stochastic"  filtered  Poisson  process.  This  three 
level  nature  of  the  signals  is  illustrated  by  Figure  1. 

1.2.2  Classical  Detection  Methods  - Gaussian  Signal  Regime 

Presently  accepted  burst-counter  and  frequency- tracker 
LV  processors  were  developed  by  analogy  with  wide-band  fre- 
quency modulation  (FM)  and  radar  receivers.  For  FM  and 
radar  applications  frequency  detection  (zero-crossing)  cir- 
cuits generally  require  about  10  db  signal  power  to  noise 
power  ratio  within  the  bandwidth  of  the  system  filters.  The 
signals  in  such  cases  are  continuous  and  the  noise  is  sta- 
tionary and  Gaussian.  Comments  by  several  speakers  at  the 
1975  Minnesota  LDA  Symposium  indicated  that  they  had  experi- 
mentally determined  that  their  burst  counter  LV  Bignal 
processors  failed  significantly  when  the  signal  power  to  noise 
power  ratio  (during  a burst)  was  less  than  10  db.  The  theory 
in  this  report  shows  that  this  10  db  condition  occurs  when 
there  are  approximately  10  photoelectron  pulses  per  elec- 
tronic response  time  x^.  This  response  time  is  rise  time  or 
pulse  width  in  the  case  of  a low-pass  filter.  We  have  also 
defined  this  signal  level  as  the  lowest  value  of  the  "Gaussian" 
signal  regime  wherein  the  photomultiplier  current  can  be 
modeled  as  the  classical  signal  plus  nonstationary  Gaussian 
noise?  (although  neither  Snyder  [6]  or  Papoulis  [7]  give 
any  helpful  rules  as  to  when  this  asymptotic  approximation 
is  valid) . For  lower  signal  power  the  signals  mUBt  be  treated 
as  Poisson. 

1.2.3  Photon  Correlation  Techniques  - Very  Low  Level  Signals 

A radically  different  approach  to  LV  signal  detection 
was  take--  by  Pike,  Oliver,  Jakeman , and  others.  Photon 
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counting  techniques  were  developed  for  use  with  low-level 
photon  resolved  signals.  The  summary,  results  of  several 
years  of  development  of  the  single-clipping  real-time 
Malvern  photon  correlator  were  described  by  Oliver  and 
Jakeman  in  a recent  book  [8] . Dr.  Pike  described  the  appli- 
cation of  photon-correlation  to  the  processing  of  LV  signals 
at  the  1972  Purdue  conference  [9] . The  presentation  was 
apparently  not  received  well  by  many  American  attendees  and 
little  has  been  done  in  the  United  States  with  the  development 
of  photon  counting  techniques  until  recently.  Increased  inter- 
est was  shown  by  attendees  of  the  1974  Purdue  conference. 

One  reason  that  the  Malvern  correlator  has  been  slow 
to  receive  acceptance  in  this  country  iB  that  the  original 
theory  for  its  use  was  based  on  the  assumption  of  many  scatters 
in  the  probe  volume  with  the  central  limit  theorem  invoked 
to  render  the  statistics  of  the  scattered  electric  fields 
Gaussian.  Since  this  assumption  is  known  not  to  predominate 
in  many  applications  of  fringe-type  LV's  in  unseeded  air, 
none  of  the  quoted  theory  for  a single-clipping  correlator  is 
directly  applicable.  Other  problems  with  the  Malvern  corre- 
lator for  high-speed  air  flow  besides  its  speed  (minimum 
time  resolution  of  50  nsec)  is  the  lack  of  any  straight- 
forward ways  to  extend  the  concept  to  the  determination  of 
flow  time  statistics  such  as  correlations  and  power  spectra. 

1.3  SCOPE 

In  Section  2.0  of  this  report  we  introduce  a 
three  level  Poisson  signal  model  for  dual  scatter  signals, 
define  signal  regimes,  and  provide  instantaneous  signal  and 
noise  equations  valid  for  al±  signal  levels.  Section  3.0 
describes  the  basic  theory  for  photon  correlation  of  LV 
signals,  and  defines  and  analyzes  a photon  counting  fre- 
quency discrimination  technique  which  we  have  called  "Dual 
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correlate  and  subtract  mode."  Section  4.0  describes  the 
design  and  the  functional  operation  of  a new  high-speed 
experimental  hardware  system.  Experiments  which  we  have 
performed  in  the  characterization  of  PMT's,  in  the  meas- 
urement of  burst  amplitude  distributions,  and  in  the  demon- 
stration of  feasibility  of  the  new  processor  are  all  described 
in  section  5.0.  Sections  6.0  and  7.0  are  devoted  to  Discus- 
sion and  Conclusions;  Appendix  A provides  a background 
information  on  Poisson  processes.  Appendix  B is  a brief 
report  in  itself  which  relates  the  burst  amplitude/rate  dis- 
tributions to  their  optical  sources  and  to  the  effects  in 
photon  counting  and  burst  counter  processors. 
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2.0  STATISTICAL  THEORY  OF  DUAL  SCATTER  SIGNALS 

2.1  POISSON  MODELS 

2.1.1  The  Signal  Current 

The  signal  current  from  a photomultiplier  tube  (PMT) 
is  modeled  as  inhomogeneous  filtered  Poisson  random  process 
(see  Appendix  A and  reference  11)  given  by 

CD 

i < t ) - l egih(t-ti)  (2.1) 

— oo 

where  t^  ■ random  time  of  occurrence  of  ith  photoelectron 
e ■ electronic  charge 

g * random  charge  gain  of  PMT  multiplier 
h(t)  ■ impulse  response  of  PMT/filter  system 
The  system  response  h(t)  is  obtained  as  a convolution  inte- 
gral of  the  PMT  impulse  response  hp(t),  the  transmission 
line  impulse  response  ht(t),  and  the  linear  filter  impulse 
response  hf(t) 


h(t)  * hp (t) *hfc (t) "hf (t) 


(2.2) 


where 


CD 

f(t)*g(t)  = / f (a)g(t  - ot)da  (2.3) 

— 00 

The  superposition  assumes  operation  in  the  linear  range  of 
the  PMT  electron  multiplier.  Th*'  use  of  the  function  hp(t) 
assumes  that  all  single  photo-electron  pulses  have  the  same 
shape  except  for  amplitude.  This  neglects  minor  random  shape 
variation. 
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The  quantity  which  relates  i(t)  to  the  classical 
optical  power  is  the  statist  cal  mean  rate  (t)  of  occur- 
rence of  the  randomly  occurring  photoelectron  pulses.  Thus 


X (t) 


HP  (t) 
hv 


(2.4) 


where 

n » product  of  cathode  quantum  efficiency  and  the 
dynode  collection  efficiency 

hv  * Photon  energy 

P(t)  « Classical  optical  power,  including  background 

light  and  a constant  component  for  dark  current. 

The  effects  of  dark  current  are  included  by  adding  an 
equivalent  power  P^.  The  model  could  be  made  more  exact  by 
adding  a separate  dark  pulse  summation  with  a separate  dis- 
tribution of  amplitudes  gd^  which  are  distributed  somewhat 
differently  than  g^;  but  this  distinction  will  not  often  be 
required. 


2.1.2  Superposition  of  Classical  Single  Burst  Signals 


The  previous  material  includes  little  which  restricts 
it  to  LV  signals.  We  now  consider  the  form  of  A (t)  which  1b 
also  treated  as  a filtered  Poisson  process. 

Rigorous  electromagnetic  theory  analysis  of  the  scat- 
tered fields  from  more  than  one  scatterer  in  the  probe  vol- 
ume shows  [12]  mixing  terms  in  P(t),  the  classical  power 
incident  on  the  PMT.  However,  in  typical  dual-scatter  sys- 
tems, the  diffraction  limited  spot  size  of  the  collecting 
lens  is  much  smaller  than  the  probe  volume;  conservation  of 
energy  arguments  show  that  in  such  cases  the  number  of 
scatterers  in  the  probe  volume  may  be  much  greater  than 
unity  with  statistically  negligible  coherent  mixing,  regardless 
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of  the  quality  of  the  collecting  lens  [l] . This  is  sig- 
nificant even  for  LV  systems  which  only  trigger  on  isolated 
large  signal  bursts  because  we  must  also  include  in  the 
model  the  effects  of  smaller  scatterers  which  may  exist  at 
higher  number  density.  We  will  take  the  position  that  at 
the  PMT  the  classical  power  P(t)  is  the  superposition  of 
the  background  light  power  and  the  power  from  individual 
scatterers  without  coherent  mixing  cross  terms.  This  will 
be  acceptable  so  long  as  the  average  number  of  scatterers  in 
one  difrraction  limited  resolution  cell  of  the  receiver  is 
less  than  unity,  and  there  are  no  strong  coherent  reflections 
from  windows  or  other  optics. 

A second  consideration  concerns  the  background  light. 
Even  when  we  neglect  coherent  mixing  of  signals,  there  are 
fluctuations  in  the  classical  background  power.  Bertolotti 
[13j  provides  a good  review  of  these  effects.  Broadband 
background  sources  can  be  largely  suppressed  by  the  use  of 
narrowband  spatial  and  wave-length  filters,  but  not  always 
adequately  enough  for  measurements  from  small  scatterers. 

In  such  cases  there  are  two  considerations.  If  the  back- 
ground is  modulated  (for  example  fluorescent  lights)  the 
mean  value  Bignal  is  easily  removed  by  electronic  filters, 
but  the  non-stationary  noise  is  not.  When  the  broadband 
background  is  "steady"  there  are  actually  significant 
classical  fluctuations  at  rates  up  to  the  optical  bandwidth. 
However,  Bertolotti  shows  that  when  the  optical  filter  band- 
width is  much  greater  than  the  PMT  electronic  bandwidth,  the 
photoelectron  statistics  behave  as  though  the  classical 
fluctuations  did  not  exist  (they  are  averaged  out) . Laser 
light  scattered  from  windows  is  not  broadband  and  may  exhibit 
undesirable  fluctuations.  This  background  should  be  mini- 
mized . 
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2.1.3  The  Random  Process  1 (t ) 


With  cognizance  of  the  preceding  discussion,  we  model 
P(t)  as  the  summation  of  a constant  background  which 
includes  broadband,  laser,  and  dark  current  sources  end  an 
inhomogeneous  filtered  Poisson  signal  process. 


h\>  > (t) 

n 


p (t) 


l X . f ( t - T . , 

©O  J J 


vjr  j 


hv 

n 


(2.5) 


where 


T . 

D 

h 


occurrence  time  of  jth  scatterer  reaching  r ^ (2.6) 


random  peak  amplitude  parameter 
vector  velocity  of  the  jth  scatterer 
location  of  nearest  approach  of  the  scatterer 
trajectory  to  the  center  of  the  probe  volume 
f (t,V,r)  = normalized  optical  system  response  function 
The  notation  in  equation  (2.5)  explicitly  showB  that  in 
general  the  shape  of  a burst  (including  signal  period,  sig- 
nal envelope,  and  pedestal)  is  a function  of  the  scatterer 
vector  velocity  and  trajectory  location.  The  response  also 
has  random  amplitude  which  depends  on  both  trajectory  and 
particle  scattering  cross  section.  In  general  the  set  of 
instants  are  independent  Poisson  random  events  whose 
instantaneous  rate,  R(t),  is  statistically  correlated  with 
the  velocity  vector. 


2.1.4  Discussion  of  Model  for  the  Classical  Signal  A(t) 

Equation  (2.6)  is  cast  in  a general  form  which  obscures 
certain  details.  First  it  assumes  that  the  velocity  of  a 
scatterer  remains  constant  while  in  the  probe  volume  with  a 
value  V(Tj,r^).  The  extended  theory  of  filtered  Poisson 
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processes  is  sufficiently  general  to  encompass  the  fact 
that  the  functional  form  of  the  optical  response  function 
f(t,V,r)  depends  on  two  vector  random  variables.*  However, 
Snyder  f6]  assumes  that  the  vector  random  parameters  are 
independent.  We  are  not  certain  at  the  present  time  what 
the  statistical  dependence  of  the  rate  function  R(t)  on  the 
velocity  V(t)  may  imply.,  We  anticipate  that  in  the  future 
other  work  related  to  Snyder's  book  will  be  applied  to  equa- 
tion (2.6)  with  less  and  less  restrictive  assumptions  con- 
cerning the  parameters. 

2.2  CONDITIONAL  SIGNAL  STATISTICS 


At  times  the  models  for  the  systems  analysis  problem 
may  be  simplified  until  analytical  methods  are  applicable. 

In  these  cases  the  use  of  conditional  statistics  will  usually 
simplify  the  analysis.  Papoulis  [7]  discusses  the  use  of 
conditional  statistics  at  length.  We  illustrate  this  tech- 
nique in  an  application  of  the  theory  in  a later  section. 
Basically  for  a multilevel  random  process  the  technique  con- 
sists of  assuming  the  higher  level  random  processes  are 
known  and  deterministic,  evaluating  conditional  expectations 
assuming  the  higher  level  processes,  then  evaluating  the 
expectation  of  the  result  with  respect  to  the  higher  level 
processes.  First  we  will  consider  statistics  of  ,i(t)  assum- 
ing the  classical  optical  signal  X(t)  is  known. 

2.2.1  Instantaneous  Mean,  Variance,  Autocovariance 

The  results  in  Appendix  A ma^  applied  to  determine 
the  instantaneous  mean,  v iance  and  ...  o-covariance  of  the 


* 

Elementary  shot  noise  theory  is  restricted  to  an 
impulse  response  function  which  is  constant  in  shape. 
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signal.  These  are  given  in  terms  of  the  function  A (t) . 

The  results  are  as  follow: 

rrn(t)  = <i(t)>  = e<g^> A (t) *h (t)  (2.7) 

a*  (t)  = < (i  (t)>-  <i(t)>)2>  = e?<gi2>A  (t)  *h2  (t)  (2.8) 

Cii^tl't2^!=  <i  (tj.)  i (t2)  > ~ <i  (t^)  ><i  (t2)  > (2.9) 

00 

= c2<g^2>/  A (a) h (t^  - a) h (t2  - a) da 
00 


where  < > denotes  statistical  expectation  and  where  the 
asterisk  again  denotes  the  convolution  integral.  These 
results  include  the  specification  that  h (t)  , the  impulse 
response  of  the  PMT  anode,  have  unit  weight,  i.e., 


oo 

/ h (t)dt  * I (2.10) 

-oo  P 


in  order  to  maintain  conservation  of  charge.  The  functions 
hfc(t)  and  h^(t)  may  include  amplification  or  loss  factors 
and  need  not  have  unity  weight. 


2.2.2  Conditional  Noise  and  SNR 

The  concept  of  signal-to-noise  ratio  arose  in  communi- 
cations theory  when  the  "noise"  was  an  additive  stationary 
Gaussian  random  process  totally  characterized  by  a mean, 
mean-square  deviation  (variance) , and  a power  spectral  den- 
sity. The  ratio  of  the  peak  or  average  signal  to  the  rms 
noise  was  a useful  measure.  The  preceding  equations  show 
the  mean-square  deviation  (variance)  to  be  an  instantaneous 
time  function  which  is  related  to  the  classical  signal. 
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Observation  of  real  signals  or  computer  simulations  such  as 
those  shown  in  Figure  2 show  that  the  concept  of  signal-to- 
noise-ratio  is  not  an  adequate  figure  of  merit  without 
careful  specification. 

As  an  example  SNR  definition,  we  consider  a low-pass  PMT 
impulse  response  as  a rectangular  function: 


h(t) 


— Rect  (t/t.) 
h n 


where 


(2.11) 


Rect  (t)  = 1,  | t | < 0.5 


(2.12) 


= 0, 


t > 0.5 


If  we  now  also  assume  that  * < < T where  T is  the  signal 

n o o J 

period  of  interest,  then  we  could  obtain  an  instantaneous 
SNR  from  equation  2.7  as 

SNR ( t)  = m.2  <t)/a.2  <c)  = <g.>V(t)/<gi2>X(t)/Th  (2.13', 


<g  > 


^i  > 


17  x(t)Th 


2 2 * 
where  the  quantity  <g^>  /<g^  > is  typically  between  0.5  and 

1.0  with  magnitude  depending  on  the  relative  variance  of  the 

PMT  single  photoelectron  pulse  gain.  For  an  ideal  tube  the 

quantity  X(t)rh  would  be  the  instantaneous  SNR.  This  is  not 

so  useful  since  it’s  a time  function  instead  of  a number. 

As  an  alternative,  wti  may  take  the  local  time  average  of 

the  SNR  given  by  (2.13)  over  a signal  cycle  near  the  peak  of 

the  pedastal.  This  would  give 
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SNRAVpeak  * Vh  (2-14> 

where  X.,  is  the  peak  value  of  the  pedastal  of  the  jth 

j 

signal  burst,  if  we  assume  sparse  non-overlapping  bursts. 

We  observe  that  equations  (2.7)  through  (2.9)  are  valid 
when  h(t)  is  a bandpass  function,  but  (2.14)  is  then  mean- 
ingless unless  we  redefine  for  a bandpass  h(t).  Also  we 
note  that  this  definition  would  be  necessary  for  meaningful 
use  with  a burst-counter  processor,  since  it  is  the  bandpass 
filtered  AC  signal  to  wide-band  noise  power  that  is  signifi- 
cant in  chat  case.  We  leave  the  subject  of  SNR  with  the 
advice  that  one  obtain  a careful  definition  before  comparing 
notes  with  others. 

2.2.3  Signal  Regimes 

The  idealized  quantity  SNRftv  peflk  given  in  equation 
(2.14)  as  XjT^  is  at  least  a useful  quantity  in  defining  a 
classification  system  of  signal  regimes  for  a low-pasB  filtered 
signal.  The  following  definitions  of  a photon-resolved,  a 
photon-limited,  and  a Gaussian  signal  regime  have  been  some- 
what arbitrarily*  identified: 

The  signal  is  photon  resolved  if  Xt^  <<  1.  In  this  case 
the  probability  of  two  or  more  photoelectron  events  occurring 
within  the  response  time  is  small.  I to  appearance  is  that 
of  individual  pulses  which  vary  in  height  due  to  the  random- 
ness of  g^.  Photon  counting  methods  are  appropriate.  The 
conditional  mean  value  of  i(t)  is  still  proportional  to  X (t) . 
ThiB  condition  is  illustrated  by  the  lower  trace  in  Figure  2. 

For  Xi^  >>  1,  the  signal  i(t)  is  asymptotically  a 


* 

See  Papoulis  [7j  page  571.  No  specific  limit  on  the 
magnitude  of  Xt^  is  given.  We  define  the  photon  limited  regime 

as  0.1  Xt.  < 10 . 
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non-atationary  Gausaian  Process.*  In  this  case  the  first 
and  higher  order  probability  density  functions  for  i(t)  at 
any  set  of  instants  (t^/t^,...)  may  be  determined  immediately 
by  plugging  the  mean,  variance  and  auto-covariance  *rom  the 
preceding  equations  into  well  known  Gaussian  formulas.  Under 
these  same  conditions  the  signal  display  appears  to  the  eye 
as  a classical  signal  m^(t)  plus  Gaussian  noise.  This  con- 
dition is  illustrated  by  the  upper  trace  in  Figure  2.  The 
major  difference  between  this  case  and  that  of  classical 
communications  theory  problems  is  that  the  a value  for  the 
noise  is  signal  (time)  and  system  dependent.  Usually,  sig- 
nals in  the  Gaussian  regime  are  suitable  for  processing  by 
classical  methods  (burst  counter  and/or  tracker) . 

The  photon  limited  regime  is  that  for  which  Xt^  is 
within  an  order  of  magnitude  of  unity.  No  mathematical 
simplif icaticns  are  possible.  Visually  the  signal  appears 
as  shown  in  the  middle  trace  of  Figure  2.  The  upper  limits 
of  photon  counting  techniques  and  the  lower  limits  of  con- 
ventional techniques  both  fall  in  this  range. 

2.3  UNCONDITIONAL  STATISTICS 

2.3.1  Long  Time  Mean,  Variance,  Autocovariance 

Equations  2.7  - 2.9  include  the  assumed  deterministic 
classical  signal  X(t)  which  is  proportional  to  instantaneous 
optical  power.  When  we  wish  later  to  evaluate,  the  long-time 
average  result  which  accumulates  during  a photon  counting 


See  Papoulis  [7]  page  571. 
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experiment,  it  will  be  necessary  to  treat  X (t)  as  an  ergodic 
random  process  with  long-time-average  equal  to  the  uncondi- 
tioned statistical  mean: 


<X  (t)  > - X 

We  also  make  use  of  the  autocorrelation  of  X (t) : 

<X (t) X (t + t) > - Rx (t)  (2.15) 

Now  from  equation  2.7  taking  the  expectation  with  respect 
to  X(t)  gives  the  mean  current  as 

00 

<i>  « e<g.><X>  / h(t)dt  (2.16) 

^ -oo 

where  the  integral  is  unity  unless  h(t)  includes  preamplifi- 
cation or  attenuation  external  to  the  PMT . In  order  to 

2 

determine  the  long  time  variance  o.  we  do  not  get  the  cor- 

1 2 

rect  answer  by  taking  the  expectation  of  cr  (t)  given  by 

equation  2.8.  Rather,  one  determines  the  conditional  value 
2 

of  <i  (t)->  by  adding  the  square  of  2.7  to  2.8.  The  expec- 
tation with  respect  to  X (t)  follows;  finally,  the  square  of 

equation  2.16  is  subtracted  from  the  unconditional  expecta- 
2 

tion  of  i (t) . When  all  these  steps  are  completed,  and 
similar  ones  for  the  unconditional  autocovariance,  we  obtain: 

OO  00 

a2  » e2[<gi>2  / Rx  (ci)f  (o)dtt  + <gi2>/  h2  (a)  da  (2.17) 

00  — 00 

00 

- <gi>2<*>2  ( / h (a) da) 2 ] 

— 00 

Cii(T)  = e^t<gi>2Rx(T)*th(i)  + <gi2><X>fh(x)  (2.18) 

oo  " 

- <g  >2<X>2(  / h (a) da ) 2 ] 

00 

2^ 
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where 


fv.(T)  = / h(a)h(a+T)da 


<2.19) 


The  second  term  in  the  expression  for  C^(t)  vanishes  for 
t greater  than  the  impulse  response  time  for  the  PMT  and 
filter  combination;  the  last  term  is  the  square  of  the  mean; 
the  first  term  is  the  correlation  of  X (t)  smoothed  by  the 
correlation  of  h(t)  with  itself. 


2.3.2  Ideal  Photon  Correlation 

An  idealized  photon  correlator  counts  all  photoelectron 
emission  events  during  successive  uniformly  spaced  clock 
periods  of  duration  At.  The  number  sequence  {n^)  which 
results  is  algebraically  manipulated  to  yield  the  summation 
of  terms  nknk+p-  In  evaluating  the  expected  value  of  the 
result  of  accumulating  such  a sum  we  encounter  the  need  to 
evaluate  the  quantities  <n^>  and  <I\+njt+p>*  These  expecta- 
tions may  be  evaluated  using  equations  2.16  - 2.19  by  assuming 


h(t)  - Rect(t/AT) 


(2.20) 


g,  - 1/e 


where  Rect(t)  was  defined  in  equation  2.12  and 


/ h (t) dt  - At 


(2.21) 


With  these  assumptions,  i(t)  is  equal  to  the  number  of 
photo-electron  events  in  the  interval  (t  - At/2,  t+Ax/2) 
and  the  formulas  reduce  to 


<i>  - <nk>  - <X>At 


(2.22) 


•**¥*‘1 
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= Oj2  = <n2>-<nk>2 * At  / Rx  (a)  A (^-) da  (2.23) 


2 2 

+ <X>At-<X>  At 


<nkW  " cn(pAT)  + <L> 


(2.24) 


ATRx(T)*A(T/AT)|T-pATfp,0 


where  the  correlation  integral  of  equation  2.19  produces  a 
triangular  function,  i.e. 


/ Rect  (—)  Rect(~)da  - AtA(t/At)  (2.25) 


where 


A(t/At)  » 1 -|t|/At  | t | <_  At 


(2.26) 


- 0 | t | > At 

When  At  is  much  smaller  than  a characteristic  signal 
period,  then  the  last  two  formulas  simplify: 


<nknk+p>  “ At  Rx(pAx)'  P ^ 0 
2 

■ At  R, (0)  + <X>At,  p ■ 0 


(2.27) 


i.e. , 


var  r^-o^  - <XAt>  + At2  (<X2>  - <X>2)  (2.28) 


var  n^  - <njc>  + At  var  ^ 


t 
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We  observe  that  the  generally  accepted  result  that  the  pho- 
ton correlation  is  shaped  like  the  correlation  of  the  classi- 
cal signal  is  true  subject  to  the  jump  discontinuity  at  zero 
delay  and  subject  to  the  triangular  weighting  function  which 
behaves  as  a low-pass  filter  with  respect  to  the  details  of 
R,  (t). 

2.4  AMPLITUDE/RATF.  STATISTICS  FOR  CLASSICAL  BURST  SIGNALS 

The  statistical  distribution  of  the  classical  burBt 
amplitudes  and  the  rate  of  occurrence  versus  amplitude  are 
very  significant  in  the  characterization  of  any  LV  signal 
processor.  It  is  generally  accepted,  for  example,  that  the 
optimum  rate  of  occurrence  of  bursts  for  a burst  counter  is 
less  than  the  inverse  burst  duration  (non-overlapping 
bursts.)  It  is  also  generally  known  that  the  error  check 
circuits  cause  a burst  counter  to  emphasize  larger  ampli- 
tude (good  signal  to  noise  ratio)  signals.  On  the  other 
hand,  a photon  counting  processor  must  emphasize  the  lower 
amplitude  signals  in  a distribution;  the  higher  amplitude 
signals  would  produce  only  a sing]e  threshold  crossing  and 
otherwise  be  neglected  by  the  system.  It  is  therefore  not 
possible  to  compare  two  different  types  of  signal  processors 
without  knowing  the  signal  amplitude  distributions  and  the 
processor  behavior  aB  a function  of  burst  amplitude  and  other 
factors.  Finally,  in  order  to  relate  processor  behavior  to 
a specified  particle  size  distribution,  one  must  first  relate 
the  particle  size  distribution  to  the  burst  amplitude  dis- 
tributions and  then  do  all  the  other  things  already  dis- 
cussed. 

During  the  initial  phases  of  this  contract  we  formu- 
lated analytical  models  and  defined  amplitude/rate  functions. 
This  task  proved  to  be  more  complicated  than  we  would  have 
imagined.  The  results  we  obtained  during  that  portion  of 
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the  study  are  not  general  in  nature  and  need  additional 
refinement  for  improved  conciseness  and  clarity.  However, 

! the  subject  matter  is  significant  for  future  LV  application 

| improvements,  so  we  have  included  the  material  from  our 

< earlier  interim  report  as  Appendix  B. 

I 
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3.0  CONCEPTUAL  FEASIBILITY  OF  A HEW  LV  SIGNAL  PROCESSOR 


In  this  section  we  postulate  a new  type  of  photon-count- 
ing LV  signal  processor  and  present  a feasibility  analysis 
which  describes  the  basic  principles.  We  begin  by  using  the 
theory  of  Section  2.0  to  derive  the  expected  result  for  an 
ideal  photon  correlator.  This  theory  is  then  extended  to 
the  Dual-Correlate-and-Subtract-Technique  which  we  have 
proposed . 

3.1  SIGNAL  MODELS  FOR  R^t) 

The  form  of  equation  2.5  is  quite  general.  In  most 
practical  systems,  the  function  f(t,Vj,tj)  which  describes 
che  optical  response  with  respect  to  particle  position  and 
velocity  is  complicated  when  the  effects  of  limiting  pinhole 
apertures  are  included.  For  the  present  we  assume  a simpli- 
fied low  turbulence  model  which  assumes  a burst 
with  perfect  contrast  and  constant  shape: 


A (t)  = A,  + ) A . f (t  - T . ) [1  + COS  ai.  (t-T,)J  (3.1) 

b L ] 1 3 J 


where 

Aj  = random  burst  podastal  amplitude  (3 

A.  = background  rate 

f (t)  - low-pass  waveform  with  peak  equal  unity 
Tj  = occurrence  time  for  burst 

ui.  = radian  frequency  proportional  to  one  velocity 
^ component  of  the  particle. 

We  may  write  equation  (3.1)  as  the  sum  of  a constant, 

low-pass  process,  A (t),  (the  pedastals) , and  a bandpass 

x,p 

process  A^ft).  Then 
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X(t) 


xb  * xlp(t) 


+ X„  (t) 
bp 


We  assume  that  there  are  several  fringes  in  the  probe  vol- 
ume so  that  the  spectra  of  X^p(t)  and  Xfap(t)  are  non-over- 
lapping. Now  Xbp(t)  is  a zero-mean  process,  and  X^p(t)  and 

X.  (t)  are  uncorrelated.  Thus  we  obtain 
bp 

VT)  ■ ♦ cXip(T)  + CXbp(T)  (3-3) 


We  now  use  the  low  turbulence  assumption  and  further  assume 
that  the  scatterers  are  uniformly  dispersed  in  space  so  that, 
R (t ) , the  rate  of  burst  arrivals  is  a constant  R.  The 
results  of  Appendix  A can  be  applied  to  derive  expressions 
for  the  three  terms  in  equation  (3.3),  the  result 


<X  (t) > ■ 


cxip<*> 


cXbp(Tl  ■ H<x 


X.  + R <x.>/  f ( t)  dt 
® J 

R<X?>/  f (t)  f (t  + t)dt 

J -fO 

go 

i <COS  U)T>/  f (t)  f (t+T  ) dt 

J ~ *•  oo 


(3.4) 

(3.5) 

(3.6) 


where  the  expectation  of  the  <cos  idt>  term  is  with  respect 
to  the  random  variable  . The  derivation  requires  that  we 
expand  the  product  of  cosines  with  the  sum  and  difference 
formula  and  approximate  the  integral  of  the  product  of  a 
low-pass  term  and  a bandpass  term  as  zero. 

We  now  assume  for  this  derivation  that  the  turbulence 
is  small  and  Gaussian  with  mean  radian  frequency  tom  and  rms 
deviation  o . Then  by  direct  application  of  the  definition 
of  expectation  we  obtain* 


In  general  <cos  gjt>  is  simply  related  to  the  statis- 
tical characteristic  function  for  the  random  variable  u>. 
since  this  function  is  defined  as  a Fourier  transform  of  the 
probability  density  function  [73. 
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co  -(ai-w  )2/2o2 

'cos  oot>  - / — e 

- <*>  /2  IT  0 

to 

-oV/2 

» e w cos  u)  t 
m 


cos  u)i  dm  (3.7) 


A simplified  expression  for  the  autocorrelation  of  the 
classical  signal  is  thus 

Rx(i)  “ {X  + R<\,>  j f(t)dt)  (3.8) 

2 2., 

->  , -0  T /2  « 

+ R<  if  > ( 1 + 4 e “ cos  0)  t)  / f ( t)  f (t+T  ) dt 

* — cn 


This  result  shows  that  regardless  of  the  shape  of  the  enve- 
lope function  f(t),  the  autocorrelation  function  has  a cosinu- 
soidal variation  at  the  mean  signal  frequency  providing  all 
the  previous  assumptions  are  met. 

The  approach  taken  by  Pike  and  others  has  been  to  fur- 
ther assume  that  the  envelope  function  f (t)  i3  also  Gaussian 
to  complete  the  derivation  of  the  shape  of  the  correlation 
function.  In  this  case  we  assume 


f (t) 


. 2 , 2 

-t  /a 


(3.9) 


where  a is  the  1/e  half  width  of  the  envelope  and  obtain 


ao 

/ f ( t)  dt  = /tt  a 

•*  on 


(fi  2 2 
/ f (t)f  (t+Odt  = ue'T  /2a 

ao 

The  final  simplified  expressions  for  the  fir 
of  X (t)  are 


(3.10) 


M, 


st  two  moments 
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<\>  = <X ( t) >=  X.  + R<X.>/jT  a 

b j 


(3.11) 


Rx(t) 


<X>2  + R<)j> (1  + j e 


-a2 t2/2  .fir 

U ' COB  umx>^2a 


2..  2 

-T  /2a 


(3.12) 


where 


Xfa  ■ background  count  rate 
R a scatterer  arrival  rate  \ 

Xj  * random  pedastal  height  from  j scatterer 
a * 1/e  half  burst  duration 

o ■ rms  deviation  of  radian  frequency  due  Vo  turbulence 


ojm  “ mean  radian  frequency  (o^/w^  <<  1)  \ 

t « delay  variable  of  autocorrelation.  \ 

The  shapes  of  typical  correlation  functions  foA  zero 
turbulence  and  10%  turbulence  are  illustrated  in  Figure  3. 
We  observe  that  functions  of  approximately  this  shape  have 
been  derived  by  others.  We  note  that  our  derivation  Bhows 
exactly  where  approximations  and  assumptions  are  U8ed,and 
also  reveals  errors  in  the  results  by  Pike  et  al.  [l5j,who 
omitted  the  Gaussian  term  resulting  from  the  pedastal. 


3.2  IDEALIZED  PHOTON  CORRELATION  OF  LV  SIGNALS 

The  number  is  the  number  of  photoelectron  emission 
eventB  in  the  interval  which  extends  ±At/2  from  the  instant 
At.  hn  idealized  photon  correlator  produces  and  sums 
delayed  products  from  the  uniformly  spaced  sequence  (n^). 

We  assume  the  total  number  of  products  accumulated  is  N. 

The  accumulator  produces  a sum  at  the  delay  value  pAt 
given  by 


" Ji  Vfc+P 


(3.14) 
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The  unconditional  expected  value  of  this  sum  is  obtained  by 
interchanging  expectation  and  summation  as 


N<n.  n.  , > 

k k+p 


(3.15) 


- NAT  (pAT) , p / 0 
= N [ <A>At  + At^R^  (0)]  , p « 0 

where  the  last  two  lines  of  equation  3.15  are  obtained  from 

equation  2.27  under  the  assumption  that  At  is  small  compared 

with  che  details  o':  R^(x).  We  see  that  the  zero  delay  term 

differs  from  the  mean  square  value  of  X by  a predictable 

amount.  The  shape  of  the  expected  result  is  otherwise  given 

by  the  function  F.^(x)  which  may  be  interpreted  via  models 

such  as  that  given  in  equation  3.12.  When  At  is  not  small 

compared  with  the  details  of  R^(t),  then  the  result  is 

smoothed  according  to  equation  2.24.  The  difference  between 
2 

the  At  multiplier  in  2.27  and  the  At  multiplier  in  2.24 
occurs  because  the  area  under  tne  triangular  function  A(t/At) 
is  just  At. 


3.3  A PHOTON  COUNTING  FREQUENCY  DISCRIMINATOR 


3.3.1  Basic  Concepts 


We  will  let  the  quantity  M be  defined  as  follows 

pg 


M = 7 m, 

kii  11 


(3.16) 


'(X^-pHiiVk-q) 


where 


-A- 


A E DC-T  R-76-81 


in,  - n,  (ri.  - n,  ) 
k k k-p  k-q 


(3.17) 


It  is  straight  forward  to  sjiow  that,  except  for  a few  end 
terms  which  do  not  affect  the  result,  M is  mathematically 
identical  to  the  quantity  & - where  R and  ft  are  defined 

by  oquatton  3.14.  It  is  for  this  reason  that  we  will  lable 
the  approach  as  the  "Dual  Correlate  and  Subtract"  technique 
even  though  the  delay  values  are  negative  instead  of  posi- 
tive.* It  is  our  purpose  to  now  demonstrate  that  the 
expected  value  of  behaves  as  a frequency  discriminator 
as  illustrated  in  Figure  4 under  conditions  which  we  will 
identify.  The  adjustment  of  the  system  clock  period  At  leads 
to  a null  in  the  expected  accumulator  value.  Measurement  of 
At  provides  a direct  measure  of  the  mean  signal  frequency  an 
we  shall  now  show  with  our  simplified  signal  models. 

From  equation  3.15  we  obtain  the  expected  value  of 
the  accumulator  as 

<Mpq>  * NAT2[Rx(pAr)  - R^qAx)]  (3.18) 


The  complete  expression  is  obtained  for  our  simplified  sig- 
nal model  by  using  equation  .3.12  in  its  entirety.  Here  we 
assume  that  pAx  and  qAx  are  both  small  compared  with  a so 
that  the  pedastal  terms  cancel  as  well  as  the  steady  term. 
This  leaves 


NAx‘R<X2> 


(3.19) 


1,1.  2,  ...  2 1,1,  2 , , . , 2 
" 2 (pAT)  " 2 V (c3Ax) 

e a cos  u>  pAx  - e a cos  w qAx 

m m . 


Rx(t)  = Rx(-i) 

The  negative  delay  implementation  v.as  more  suitable  for 
the  hardware  design. 
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Equation  3.19  appears  to  be  complex.  Now  we  choose  to 
require  that 


p - 3q  (3.20) 

and  observe  the  behavior  of  equation  3.19  near  the  values 
of  At  specified  by  letting  qAx  ha  one  quarter  of  the  signal 
period  where  both  cosine  terms  will  vanish. 

qAt  - 2 Tr/4aj  « T/4  (3.21) 

^ 'mm 

The  shape  of  the  term  in  braces  is  plotted  in  Figure  4 
under  the  assumption  of  many  fringes  in  the  probe  volume 
(a  large)  and  low  turbulence  (o  small).  Thus  the  quantity 
plotted  is  simply 

[cos  3qAt  uj  - cos  qAt  to  ] (3.22) 

m ’ m 

Figure  4 illustrates  the  expected  behavior  of  the  accumula- 
tor sum  for  changes  in  the  mean  signal  frequency  ui  . If 
the  system  clock  frequency  is  changed  to  change  At,  then  the 
response  is  the  product  of  At  and  the  curve  shown  in  the 
figure.  The  shape  of  the  curve  is  affected  but  the  zero 
crossing  locations  are  not. 

3.3.2  Selection  of  Delay  Constants 

The  theory  of  Section  3.3  does  nor  unique], y specify  the 
relationship  of  At  ( the  system  clock  period)  to  Tm  the  sig- 
nal period.  For  a given  value  of  signal  period,  T , the 
largest  possible  value  of  At  for  an  acceptable  null  is  when 
q = 1 and  At  - T /4.  This  produces  the  least  variability 
error  (see  section  6.0)  and  the  most  bias  error  due  to  time 
smear  (the  triangular  weighting  function  in  equation  2.24). 
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The  value  q = 1 also  allows  the  highest  signal  frequencies 
to  be  measured  for  a given  maximum  system  clock  frequency. 

The  bias  errors  should  be  reduced  at  the  expense  of 
increased  variability  error  and  reduced  maximum  signal  fre- 
quency by  using  q - 2,  At  ■ Tm/8.  In  order  to  facilitate 
experimental  research , our  design  for  a research  instrument 
allows  selection  of  p and  q over  a wide  range. 
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4.0  HARDWARE  PHOTON  PROCESSOR 


4.1  DESIGN  RATIONALE 


We  have  designed  and  constructed  a feasibility  proto- 
type photon  processor*  to  test  the  Dual  Correlate  and  Subtract 
concept  described  theoretically  in  Section  3.0.  Although 
this  was  a research  system,  it  was  designed  so  that  it  would 
outperform  the  existing  Malvern  photon  correlator  in  several 
significant  ways  regardless  of  the  feasibility  of  the  Dual 
Correlate  and  Subtract  concept.  In  this  manner,  we  have 
assured  that  the  AEDC  would  obtain  an  advanced  photon  pro- 
cessing system.  In  designing  the  system  we  have  incorporated 
the  following  guidelines: 

1.  The  design  objective  maximum  clock  speed  was  100  MHz 
(10  nsec  time  resolution)  with  a guaranteed  40  MHz 
operation . 

2.  The  system  would  utilize  a four  bit  multiplier  which 
operates  at  the  clock  frequency  (down  to  10  nsec  per 
multiply)  to  avoid  the  ambiguities  of  interpretation 
of  the  single-clipping  correlator. 

3.  The  system  was  to  be  electronically  switchable  in 
configuration  to  provide  a variety  of  different 
functions  for  different  research  and/or  measurement 
purposes,  including  direct  photon  correlator. 

4.  The  system  would  be  designed  to  be  computer  compati- 
ble for  extended  dynamic  modes. 

5.  The  system  and  development  costs  would  be  minimized 
Bubject  to  the  other  goals. 

In  order  to  satisfy  the  cost  and  speed  objectives  the 
system  was  designed  with  MECL  10000  logic  on  special  wire 
wrap  panel.  Significantly  higher  speed  would  have  required 


Detailed  design  and  construction  provided  by  GW. 
Electronics,  Atlanta,  Ga . 
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the  use  of  high-speed  printed  circuit  board  techniques  and 
higher  speed  ECL  logic.  The  use  of  printed  circuit  board 
techniques  would  have  significantly  increased  the  develop- 
ment time  and  hence  costs.  Although  TTL  circuits  and  tech- 
niques could  have  been  used  in  a straight  forward  manner  &■:, 
lower  speeds,  the  use  of  th^-se  circuits  at  100  MHz  would  have 
been  difficult  or  impossible  and  we  chose  the  more  expensive 
but  faster  MECL  10000  to  c/oid  additional  expense  in  design 
and  debugging  time. 

A synchronous  pipe  line  multiplier  design  was  imple- 
mented which  allows  the  4 bit  multiplication  to  be  obtained 
via  sequential  shifts  and  additions.  This  multiplier  tech- 
nique requires  the  use  of  binary  arithmetic.  Simplicity  and 
cost  concerns  then  dictated  that  all  high-speed  switching 
and  arithmetic  be  done  in  binary  arithmetic.  In  order  to 
minimize  cost  for  the  feasibility  study  while  max.imizing  the 
system  potential,  the  system  was  designed  so  that  it  could 
be  operated  by  a person  via  toggle  switches j and  chen,  if 
feasibility  were  demonstrated,  it  could  be  conveniently 
interfaced  to  a minicomputer  or  special  purpose  dedicated 
microcomputer-controller  for  user  oriented  ease  of  control  and 
data  acquisition.  The  potential  for  microcomputer  control 
will  allow  high-speed  reconfiguration  in  the  future  for  modes 
of  operation  not  presently  possible  such  as  rapid  sequential 
autocorrelation  and  multiple  photon  interval  statistics. 

4.2  SYSTEM  DEFINITION 
4.2.1  Overview 

The  delivered  system,  as  illustrated  in  Figure  5,  con- 
sists of  an  Ortec  9301  preamplifier,  Ortec  9302  amplifier/ 
discriminator,  Ortec  401/402m  minibin  and  power  supply,  a 
System  Donner  6151  laboratory  counter  (SDC) , a WAVETEK  1001A 
sweep  generator  (used  as  the  variable  local  oscillator) , and 
the  SAI/GW  Electronics  Photon  Processor. 
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The  SDC  counts  directly  without  prescaling  to  200 
MHz  and  has  user  installable  options  for  computer  interface 
cards  for  control  and  data  acquisition.  It  is  used  both  in 
diagnostics  and  system  control  and  data  acquisition.  The 
WAVETEK  oscillator  is  adjustable  manually  or  by  remote  pro- 
gramming analog  voltage  over  the  range  500  KHz  to  300  MHz. 

The  Ortec  preamplifier  is  required  with  a low-gain  PMT  such 
as  an  RCA  8644  and  is  helpful  with  a higher  gain  tube  such 
as  the  RCA  931A.  The  preferred  RCA  tubes  for  photon  counting 
are  the  RCA  8850  or  C31024  tubes.  The  preamplifier  is  not 
required  with  these  tubes  (Tubes  by  othar  PMT  Manufacturers 
have  not  been  investigated.  We  do  not  claim  any  superiority 
for  RCA.  RCA  tubes  were  available  during  the  experimental  work. 
When  used  it  is  significantly  important  that  the  preamp  be 
located  physically  directly  at  the  PMT.  The  combination  of 
any  of  the  above  mentioned  photomultiplier  tubes,  and  the 
Ortec  preamp  and  amp/diacriminator  haB  a maximum  speed  of 
approximately  125  MHz  with  112  MHz  being  a more  reliable 
figure.  (The  dead  time  observed  on  an  oscilloscope  with  PMT 
input  to  the  9302  unit  was  10  9 nsec;  output  pulses  were 
obtained  to  125  MHz  by  driving  the  9302  input  with  an  oscil- 
lator and  increasing  the  drive  slightly  at  about  110  MHz.) 

The  other  major  system  component  is,  of  course,  the 
procesaor  we  designed  and  constructed.  It  is  described  in 
detail  in  the  next  section/  and  illustrated  in  Figure  7. 

The  coaxial  cable  UBed  between  the  9301  and  9302  units 
and  its  terminations  are  critical  elements.  We  used  precision 
50  (2  cable  and  still  found  it  desirable  to  use  a 50£?»  2x 
attenuator  at  the  output  of  the  preamp  to  reduce  pulse  reflec- 
tions which  can  cause  false  discriminator  correlations  when 
a tube  with  large  single  photo-electron  pulse  variation  is 
used. 
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By  way  of  completeness  we  include  a potential  system 
diagram  in  Figure  6 which  illustrates  a proposed  interface 
first  with  a microcomputer  used  in  a minicomputer  configu- 
ration and  later  with  a micro  processor  controller  and  data 
acquisition  system, 

4.2.2  Input/Output 

There  are  three  BNC  input  jacks  on  the  rear  panel 
which  are  for  the  variable  oscillator  clock  input,  the  out- 
put of  the  photon  discriminator  (ORTEC  9302)  and  the  gate 
control  pulse  from  the  Rystron  Donner  Counter.  There  are 
three  output  BNC's  which  provide  decade  scaled  versions  of 
the  clock  (+  10,  > 100,  and  : 1000). 

The  front  panel  is  illustrated  by  the  photograph  in 
figure  7.  The  top  row  of  display  lights  consists  of  5 
LED's  (left  hand  side)  which  continuously  display  the  output 
of  the  4 bit  counter.  The  remainder  of  the  top  row  of  lights 
is  the  accumulator  overflow,  sign,  and  magnitude.  The  lower 
row  of  lights  display  the  states  of  the  control  address 
(4  bits)  the  control  data  word  (4  bits) , the  Read  Address 
(4  bits)  and  the  Read  data  (4  bits) . The  toggle  switches 
under  the  lights  correspond  to  the  lights  above.  The  sepa- 
rate toggle  switch  to  the  left  provide  an  "enter"  pulse  to 
enter  control  data  to  a selected  address  or  to  change  the 
read  address.  The  two  switches  at  the  lower  right  are  for 
addressing  the  control  gate  to  either  the  system  clock  or 
the  photon  data  Btream. 

4.2.3  Processor  Circuit  Description 

Details  of  the  circuits  used  in  the  photon  processor 
have  been  provided  to  the  Arnold  Center  in  the  form  of 
schematic  wiring  diagrams  which  form  a portion  of  a processor 
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Figure  6.  Proposed  Extended  System  with  Digital  Counter. 
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manual  deliveded  with  the  feasibility  processor.  Figure  8 
is  a functional  diagram  which  illustrates  all  internal 
subsystems.  The  following  paragraphs  describe  the  sub- 
system blocks  illustrated  in  Figure  8.  For  Control  and 
Read  blocks  not  discussed  below  see  the  appropriate  Control 
vocabulary  definition  in  Section  4.2.4  or  Read  Vocabulary 
definition  in  4.2.5. 

Clock  Shaping  Circuit  - The  sinusoidal  clock  input 
is  internally  terminated  with  50Q.  Shaping  is  accom- 
plished by  a flip-flop  operated  in  a monostable  mode. 

Timing  is  by  propagation  delay  and  a length  of  transmission 
line  which  is  tuned  for  optimum  clook  pulse  width. 

3- Decade  Scaler  - has  outputs  at  1/10 , 1/100,  and 
1/1000  of  the  system  clock  frequency.  Outputs  are  in  MECL 
logic  levels  and  counting  is  synchronous. 

SDC  Control  Circuit  - provides  resets  and  clears  at 
the  beginning  of  each  measurement  cycle  to  initialize  the 
machine.  The  SDC  Control  Circuit  also  causes  data  to  be 
held  while  the  TTL  input  is  low  so  that  reading  the  front 
panel  may  be  accomplished.  Two  switches  on  the  front  panel 
control  options  on  the  SDC  control  circuit.  The  switches 
are  located  at  the  lower  right  corner  of  the  front  panel 
and  are  labeled  D and  C.  When  placed  in  the  up  position  they 
inhibit  the  Data  and  Clock  respectively  during  the  time  the 
control  input  is  lov/.  In  normal  operation  these  switches  are 
left  in  the  down  position. 

Clock  Tree  - System  clock  distribution  is  by  three 
packages  of  multiple  output  gates.  This  extensive  clock 
buffering  is  necessary  in  spite  oi  MECL's  high  low-frequency 
fanout  capability  due  to  the  extreme  speed  of  the  clock  pulse 
and  the  capacitive  loading  of  the  clock. 

4- Bit  Counter  - This  is  a dual  counter,  that  is  one 
phase  counts  while  the  other  is  settled  and  connected  to 

the  comparator.  There  is  a dead  time  of  approximately  2 nsec 


Figure  8.  Subsystem  Block  Diagram  of  Photon  Processor. 
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per  clock  period  at  the  clock  transition  during  which  the 
counter  switching  occurs  and  any  data  will  be  ignored. 

The  counter  has  a programmable  stop  at  1 , 2,  4,  8,  15  and 
vill  never  overflow  into  an  erroneous  state. 

Max  Count  Control  - Like  all  control  words  in  the 
machine  Max  Count  Control  operates  by  entering  a control 
word  into  the  appropriate  address. 

Comparator  - The  Comparator  or  clipper  compares  the 
output  of  the  4-Bit  Counter  with  a number  from  the  Clip 
Leve1  Control  and  assigns  a binary  value  to  Nk  ^ based  on 
the  results  of  that  comparison.  If  the  photon  count  is 
greater  than  the  programmed  clip  level  then  N.  ^ will  be 
assigned  a value  of  1,  otherwise  Nk ' will  be  assigned  a 
value  of  C . 

Read  Xk  - causes  Xk  to  be  connected  to  the  Read  Data 
Bus  as  long  as  the  correct  address  is  held.  For  detailed 
information  on  read  addressing  and  interpretation  of  output 
data  see  faction  4.2.5. 

X Y Multiplexer  - selects  either  clipped  or  unclipped 
data  and  transmits  this  to  Yk  lines  and 
combination. 

Delay  A - This  is  a digital  delay  line  with  a length 
choBen.  by  a decoder.  All  data  passes  through  all  twelve 
delays  and  is  lost  when  clocked  out  of  the  twelfth  position. 

Delay  B - Similar  to  Delay  A,  this  delay  element  has 
taps  at  0 , 2 , 4 , and  6 . 

Subtractor  - Subtracts  the  output  of  the  Subtract  Mul- 
tiplex from  the  output  of  the  Delay  R.  Subtractor  is  on 
two's  complement  arithmetic  but  both  inputs  are  positive  so 
no  input  conversion  is  required.  Output  iB  five  bits  (4  bits 
plus  sign)  and  is  converted  to  sign  and  magnitude  represen- 
tation for  the  multiplier. 

Subtract  Multiplex  - causes  either  Yk^A  or  the  Subtract 
Constant  to  be  subtracted  from  yjc_a_b  cased  on  the  state  of 
the  Mode  Select. 


Xj^  lines  in  any 
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Subtract  Constant  - is  a register  storing  a constant 
entered  as  a control  word.  See  Section  4.2.4. 

Mode  Select  - is  two  bits  of  information  stored  in 
order  to  define  the  state  of  the  Subtract  Multiplex  and 
the  Multiply  Multiplex.  It  is  entered  as  a control  word. 

See  Section  4.2.4. 

Multiplier  - multiplies  the  outputs  of  the  Subtractor 
and  the  multiply  Multiplex.  Arithmetic  is  in  signed  magni- 
tude. The  Multiplier  is  a fully  pipelined  shift  and  add 
design. 

Multiply  Multiplex  - selects  either  or  1,  as  one 
input  to  the  Multiplier  based  on  the  state  of  the  Mode 
Select. 

Read  Accumulator  - accomplished  sequentially  by  5 
separate  Read  control  addresses  (or  visually  directly  from 
LED  display).  See  Section  4.2.5. 

Accumulator  - This  is  a pipelined  synchronous 
accumulator  in  two's  complement  arithmetic.  Conversion  of 
the  input  to  two's  complement  is  by  exclusive  or  gfites  and 
a carry  into  the  first  stage.  The  overflow  signal  is  gen- 
erated by  OR-ing  together  three  bits  more  significant  than 
the  most  significant  bit  displayed  or  available  through 
reads.  This  means  that  it  is  possible  for  inetance  to  over- 
flow positive  and  add  on  enough  negative  to  come  out  of 
overflow.  Output  conversion  to  signed  magnitude  represen- 
tation is  by  exclusive  or  gates  and  therefore  all  negative 
numbers  are  in  error  by  a one.  This  error  is  small  compared 
to  the  width  of  the  accumulator  and  its  correction  does  not 
justify  the  amount  of  circuit  complexity  which  would  be 
required . 

4.2.4  Control  Vocabulary 

Control  is  provided  by  programming  control  words  into 
various  aduresses.  The  addresses  and  definition  of  the 
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meaning  of  the  various  control  words  are  given  in  Table  1. 
4.2.5  Read  Vocabulary 

When  a Read  address  is  entered,  four  bits  of  data  are 
continuously  displayed  by  the  LED  data  lights.  For  computer 
controlled  operation,  the  clock  would  be  stopped,-  and  then 
the  various  addresses  would  be  read.  The  addresses  and  data 
definitions  are  given  in  Table  2 1 

4.3  MODEG  OF  OPERATION 

The  SAX  processor  may  be  operated  in  static  modes  at 
present.  It  has  the  potential  for  operating  in  dynamic  modes. 
A static  mode  is  one  in  which  the  machine  configuration  is 
selected  by  hand  or  by  computer  control  and  then  data  is 
taken  for  some  length  of  time.  A dynamic  mode  is  one  in 
which  the  SAI  processor  is  rapidly  re-configured  during 
the  data  collection  process  by  computer  control.  Some  of 
the  static  mcdes  arei  1.  Moan;  2.  Mean  Square;  3.  Single 
Point  Correlation;  4.  Single-Point  CroBB-Correlationj 
5.  Dual  Correlate  and  Subtract.  There  are  other  static 
configurations  and  many  possible  dynamic  modes.  Some  of 
interest  ares  sequential  autocorrelation,  counting  probability 
histogram,  diagnostic  sequences,  and  others.  We  are  estab- 
lishing here  the  procedures  and  machine  codes  and  sequences 
which  produce  the  static  modes  of  operation  just  listed. 

Our  basic  approach  for  static  modes  is  as  follows; 

1.  Establish  the  Dual  Correlate  Mode  with  a basic 
set  of  the  optional  choices. . take  data  or 
reconfigure  as  follows. 

2.  Select  other  static  modes 

3.  Select  other  optional  choices. 
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coxtroi 

ppfq-r- 

CLEMENT 

DATA 

jsml.... 

JC1.CCTION 

0000 

lU’iiERVtO  TOR  fOTURE  OPTIONS 

0001 

MIKKVSD  ron  rUTVU  OPTION! 

0010 

4-blt  counter 

0000 

•top  at  1 

0001 

Step  !t  i 

0011 

•top  ai  a 

on: 

•top  tt  0 

1111 

(tap  at  It 

0011 

Comparator 

tbod 

(abed)  » * • ellp  leval  (O-IM 
n (*)  • 1,  n > m • 0,  n S a 

A100 

XY  Multiplex 

OOkx 

X • full,  Y • full 

01  xx 

X • full.  Y • CUppud 

lOxx 

X « cUppad,  Y • full 

UN M 

X • CUppad,  Y « CUppad 

0101 

Delay  A 

0000 

Delay  • 9 aleak  period# 

0001 

Delay  • 1 alack  parlode 

0010 

Delay  • 2 eloek  periods 

oon 

Deity  • J aleck  period! 

0100 

Delay  • 4 clock  period! 

0101 

Delay  • I clock  perloda 

ouo 

Delay  » 4 aleak  period# 

1000 

Delay  • 1 alack  period! 

1001 

Daley  • 0 aleak  period! 

10\0 

Daley  • )9  cloak  period! 

1100 

Delay  • 11  aloah  period! 

QUO 

tally  • 

wcOO 

Delay  ( • 0 

xxoi 

Delay  • 1 

WJ.O 

Delay  • 4 

nil 

Delay  • 1 

0111 

Mode  (eleet 

woo 

aubtraat  eorarant,  eiultiply  by  1 

wOl 

■ubtraat  aomtant,  Multiply  by  a\ 

wlO 

aubtraet  Yk-X,  Multiply  by  1 

MX  11 

subtract  Yfc-A*  multiply  by  Xx 

1000 

aubtraction 

constant 

(abed) 

poeitive  binary  word  (O-.'J) 
equal  to  lubtraaiion  conatant 

1001  to 

mi  mmrvr9  fox  nrruxi  option* 


Table  1 Control  Vocabulary. 
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Table  2.  Read  Vocabulary 
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4.3.1  Establish  Dual  Correlate  Mode 

The  following  basic  sequence  removes  the  max  count 
limits , the  single  clipper,  and  sets  the  delays  for  eighth 
period  increments.  (Quarter  period  increments  would  be 
A *1,  B ■ 2,  instead  of  A ■ 2,  B » 4) . The  machine  com- 
putes 


% ■ nk  <nk-A-F-nk-A> 


The  code  is  entered  as  shown  in  Table  3. 


ADDRESS 

DATA 

OPERATION 

0010 

1111 

max  count  ■ 15 

0011 

0000 

clip  level  at  zero 

0100 

00XX 

full  multiply 

0101 

0010 

delay  A ■ 2 

0110 

xxlO 

delay  B - 4 

0111 

xxll 

subtract  Y^_A  and 

multiply  by 

1000 

0000 

subtract  0 

Table  3.  Dual  Correlate  Mode. 

4.3.2  Select  Other  Static  Modes 

Prom  the  basic  static  mode,  single  codes  are  used  to 
change  modes  to  the  ones  indicated  in  Table  4. 
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MODE 

FUNCTION  m. 

CONT. 

DATA 

K 

ADDRESS 

mean 

nk-A-B 

0111 

XX00 

mean  square 

nk2 

0111' 

xxOl 

0101 

0000 

0110 

XX00 

autocorrelate 

nk  nk-A-B 

set  delay 

0111 

xxOl 

(0101 

abed 

\0110 

abed 

Table  4.0  Other  Static 

Modes 

4.3.3  Clipping  and  Other  Options 

A variety  of  optional  choices  are  allowed  by  the 
processor  logic.  The  following  listing  gives  several  of 
the  most  likely  to  be  used  options. 

4. 3. 3.1  Cross  Correlate  with  Clipped  Sequence.  In  any  mode 

the  clipped  version  of  n^  (which  is  may  be  applied 

to  either  the  or  branches.  A clip  level  should  also 
bo  selected.  These  choices  are  made  with  control  addresses 
3 and  4 (0011  and  0100) . See  machine  language  Table  1. 

4. 3. 3. 2 Max  Count.  The  maximum  value  of  n^  may  bo  restricted 
to  1,  2,  4,  or  8 instead  of  15  by  using  the  desired  max  count 
number  as  the  data  word  to  address  2 (0010). 

4. 3. 3. 3 Set  Delay.  For  the  autocorrelate  mode  the  delay  is 
the  sum  of  delay  A and  delay  B.  The  data  word  for  the  A 
delay  is  equal  numerically  to  the  delay  excluding  7,  11,  13, 
14,  15).  The  data  word  for  the  B delay  is  the  log  base  2 of 
the  delay. 
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4.  3.  3. 4 Subtraction  Constant.  If  a conatant.  other  than 
zero  is  to  be  subtracted  it  should  be  loaded  into  control 
register  8 (0100) . This  feature  is  provided  to  reduce  the 
mean-square  level  stored  in  the  accumulator  when  the  auto- 
correlate  mode  is  used  and  the  steady  optical  power  compo- 
nent is  large  compared  with  the  time  fluctuating  component.) 

4.4  OPERATING  PROCEDURES 

A manual  measurement  with  the  photon  system  requires 
several  operations  with  the  Wavetek  oscillator#  the  labora- 
tory counter#  and  the  9302  discriminator.  These  are 
described  briefly  here. 

4.4.1  Discriminator  Threshold 

The  combined  gain  of  the  PMT,  the  preamp,  the  attenu- 
ation (if  any),  the  amplifier  portion  of  the  9302, and  the 
discriminator  threshold  must  be  chossn  to  satisfy  the  follow- 
ing requirements i a)  maintain  speed  b)  avoid  saturation  of 
! amplifiers  by  large  photoelectron  pulses  c)  avoid  missing 

small  photoelectror.  pulses  d)  avoid  setting  threshold  down 
into  the  amplified  Gaussian  electronic  noise  from  the  preamp 
or  the  ringing  tails  of  the  larger  single  photoelectron 
pulses.  In  practice#  these  constraints  could  not  all  be  met 
simultaneously  with  either  the  RCA  8644  or  the  RCA  931A  tube 
due  to  the  large  range  of  amplitudes  of  the  single  photo- 
electron pulses.  In  the  future  better  photon  counting  tubes 
, should  be  used. 

4.4.2  Adjustment  of  Clock  Frequency  and  Use  of  Counter 

For  an  autocorrelate  mode  the  Wavetek  oscillator  fre- 
quency is  adjuntod  by  using  the  laboratory  counter  to 
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accurately  determine  a convenient  delay  resolution  incre- 
ment. The  oscillator  is  then  left  at  that  frequency  while 
the  measurement  is  made  at  different  delays.  For  the  Dual 
Correlate  Mode,  the  clock  frequency  is  adjusted  by  an  opera- 
tor receiving  feedback  from  the  accumulator  display  lights 
until  a null  is  reached,  then  the  oscillator  frequency  is 
read  with  the  laboratory  counter  to  determine  the  clock 
frequency  accurately.  This  feedback  can  either  be  with  the 
processor  running  continuously  (best  for  very  low  photo- 
electron rates)  or  repetitively  for  fixed  duration  gate 
intervals  provided  by  the  lab  counter  output  gate.  In 
either  casts,  the  "B"  counter  of  the  lab  counter  is  counting 
its  internal  crystal  oscillator  so  that  the  "A"  counter  may 
be  used  to  count  the  total  number  of  discriminator  output 
pulses  during  the  gated  interval.  (There  are  two  isolated 
but  equal  outputs  of  the  ORTEC  9302  Amp/Disc. ) . This  simul- 
taneous count  may  be  used  for  normalization  purposes. 
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5.0  EXPERIMENTAL  MEASUREMENTS 

Four  different  sets  of  measurements  were  made  during 
the  contract  period.  The  details  of  all  of  these  measure- 
ment? have  been  documented  in  informal  monthly  progress 
reports.  Only  the  general  nature  of  the  experiments  and 
significant  data  are  summarized  briefly  here. 

5.1  PRELIMINARY  MEASUREMENTS 


This  sot  of  experiments  were  conducted  in  October 
1974  prior  to  the  construction  of  the  photon  processor. 

5.1.1  Photomultiplier  Tubes 

Several  different  types  of  observations  were  made 
with  RCA  8644,  931A,  and  G850  photomultiplier  tubes.  First, 
observation  of  single  photoelectron  pulses  on  a fast  oscil- 
loscope demonstrated  the  advantage  of  uho  8850  tube  over 
the  other  two.  When  observing  the  brightness  of  the  oscil- 
loscope display  versus  amplitude  of  the  pulses  there  is  clearly 
a preferred  amplitude  with  a ±15-25%  standard  deviation 
(qualitative  judgement)  for  the  8850  tube.  The  other  two 
tubes  produce  the  brightest  display  at  zero  amplitude, 
actually  at  the  trigger  level,  with  gradually  decreasing 

brightness  with  increasing  amplitude.  With  the  8850  tube 
it  was  easy  to  set  a threshold  trigger  level  above  the 
electronic  noise  and  bel.ow  the  minimum  pulse  amplitude  of 
most  of  the  pulses.  This  was  not  possible  with  the  other 
two  tubes.  The  pulse  amplitudes  from  the  B850  was  approxi- 
mately 100  times  that  of  the  8644  with  the  voltages  used. 
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No  preamp  was  needed  with  the  8850  tube. 

Two  RCA  8644  photomultiplier  tubes  assemblies  were 
testod  for  absolute  anode  sensitivity,  quantum/ collecting 
efficiency  product,  and  qualitative  waveform  features.  The 
first  unit  was  an  unlabled,  flat  black  assembly  which  we 
designate  Model  0.  The  second  unit  was  PMR-20,  serial  1. 

Both  units  were  tested  at  a supply  voltage  of  1800  v at 
the  green  argon  laser  wavelength  (514.5  nm)  . The  pulse 
waveform  of  the  Model  0 unit  exhibited  considerably  more 
ringing  and  therefore  leas  threshold  triggering  stability 
than  did  the  PMR-20  unit.  The  measurement  results  are  pro- 
vided below  in  Table  5.  We  note  that  while  the  anode 
sensitivity  results  are  needed  for  calibration  purposes  in 
the  burst  height  distribution  studies,  it  is  the  n quantity 
which  limits  both  photon-counting  and  burst  counter  processors 
since  variations  in  DC  dynode  gain  are  easily  compensated  by 
preamplifier  gain  and  do  not  affect  the  signal-to-noise 
ratio. 


Anode  Sens 

Amp a /watt 

RCA  typical  (spec)  16,600 
Model  zero  18,600 

PMR-20  10,600 


method  1 


(0.125) (0.8)- 
0.10 
0.035 


nF 

method  2 

0.10 

0.07 

0.04 


Table  5.  Measurement  Results 


Despite  the  clear  advantages  for  photon-counting,  the 
9850  tube  was  not  suitable  for  use  with  signals  large  enough 
for  burst-counter  operation  due  to  saturation  of  the  dynode 
divider  current.  This  could  be  remedied  by  changing  the 
dynode  divider  to  provide  more  current,  but  this  was  not 
done  during  the  contract  period.  Hence  the  poorer  tubes 
were  used  in  later  photon  counting  experiments  where  compari- 
son with  a burst  counter  processor  was  desired. 
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5.1.2  False  Triggering  of  Burst  Counter 

An  experiment  was  performed  in  which  an  ARO  Model  8 
Durst  Counter  processor  was  used  with  a laboratory  LDV 
experiment.  At  low  signal  levels  with  low  trigger  threshold 
settings  on  the  oscilloscope  used  for  triggering  the  processor, 
a significant  fraction  of  the  trigger  events  were  from 
anomalously  large  single  photo-electron  pulses  (RCA  8644 
tube  again) . Significant  improvement  in  trigger  relia- 
bility was  obtained  by  using  a low  pasB  filter  in  the  trigger 
circuit  with  bandwidth  from  2 to  4 times  the  pedastal  band- 
width. This  is  much  less  bandwidth  than  is  needed  for  the  AC  } 

signals.  The  technique  is  only  effective  when  the  mean  aig-  ] 

nal  visibility  is  high  since  the  pedastal  is  used  for  trig-  j 

gering.  Such  techniques  would  not  be  necessary  with  larger  1 

PMT  stage  gain  and  smaller  pulse  height  variation.  \ 


5.1.3  Photon  Correlation  Demonstration 


I 


I 
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Mr.  Bob  Koelzer  and  Mr.  Charles  Sheckells  of  Honeywell 
Saicor  provided  a demonstration  of  the  model  SAI  43A  digital 
correlator  with  photon  correlation  option.  A laboratory  LV 
system  was  set  up  with  an  ARO  model  8 burst  counter  processor 
for  comparison.  An  RCA  8644  PMT  was  used  in  the  experiment. 
The  experimental  arrangement  consisted  of  a forward  scatter, 
fringe  type  LDV  system  with  the.  probe  volume  located  respec- 
tively three  inches  and  one  inch  in  front  of  the  intake  of  a 
small  laboratory  wind  tunnel.  In  the  first  test  the  optical 
system  consisted  of  a small  1.6  mw  He  Ne  laser,  50-50  self 
aligning  splitter  blocks,  a 12.1  cm  fl  transmit  lens,  a pair 
of  19.2  cm  lens  in  a one-to-one  imaging  arrangement,  a folding 
mirror,  and  an  RCA  8644  PMT  with  a pinhole  and  laser  line 
filter . 
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A second  test  was  conducted  with  the  optical  system 
changed  in  the  following  manner: 

Transmit  focal  length  = 62  cm  Fringe  Spacing  = 15.4yjm 

Pinhole  diameter  - 400pm  Probe  Volume  Diameter  = 400pm 

Half-angle  between  beams  - 1.17® 

Other  system  parameters  were  unchanged  except  as  noted  below. 

In  the  conduct  of  the  second  test,  burst  counter  data 
was  first  taken  with  the  system  as  described?  then  a 10  x area 
ratio  restricting  aperature  was  placed  over  the  collecting 
optics  to  reduce  the  collecting  aperture  to  f/11.7.  Under 
these  conditions  both  burst  counter  and  photon  correlator 
data  were  taken  {12  minutes  of  burst  counter  data  followed 
by  15  minutes  of  correlator  data  (real  time  equivalent  0.375 
sec).  Data  sets  were  obtained  with  a short  (12.1  cm)  and  a 
long  (62  cm)  transmitter  configuration,  both  forward  scatter, 
using  a 1.0  mw  He  Ne  laser,  uaseeded  natural  air,  at  veloci- 
ties of  0.91  m/sec  and  5.8  m/sec.  respectively.  Burst 
counter  measurements  with  printer  output  were  also  made  for 
each  optical  geometry  and  flow  velocity.  The  mean  period  and 
rms  deviation  have  been  computed  from  60  values  of  the  burst 
counter  data  (before  lOx  attenuation)  for  comparison  with 
photographic  records  of  the  photon  correlator  display;  these 
are  shown  in  the  following  tabic.  The  data  collection  time 
for  the  correlator  iB  unknown  for  the  12.1  cm  transmit  lens 
case,  but  was  approximately  30  minutes.  The  time  for  the  62 
cm  results  was  900  seconds  (0.375  seconds,  real  time  equiva- 
lent) , with  the  good  data  rate  from  the  burst  counter  being 
0.2/sec  in  this  latter  case  (with  lOx  attenuation). 

Mean  Mean  Standard 

Transmit  Lens  Velocity  Fringe  Period  Period  Deviation 

Focal  Length P^.'lod  Photon  Counter  Counter 

12.1  cm  0.91m/sec  3.0  pm  3.2±0.1  3.29psec  0.22 

psec 

62  cm  5.8m/sec  15.4  pm  2.6±0.1  2.66psec  0.05 

r sec 

Table  6.  Photon  Correlation  Demonstration  Results. 
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The  experiment  indicated  that  the  photon  correlator 
provided  a more  sensitive  detection  technique  than  the 
ARO  processor . The  time  required  would  have  been  less  if 
tne  correlator  multiplier  could  have  been  fast  enough  for 
real  tir  e data  processing. 

5.1.4  Distributions  of  Burst  Rate  Vs.  Amplitude 

Low  speed  measurements  were  made  to  determine  the  dis- 
tribution of  LV  signal  burst  amplitudes  for  a 147  mw  back- 
scatter  Argon  LV  system  as  illustrated  in  Figure  9 . The 
parameters  of  the  system  were  chosen  similar  to  those 
expected  in  future  Arnold  transonic  windtunnel  tests. 

Absolute  quantities  were  measured  so  that  the  data  could  be 
scaled* to  typical  transonic  velocities  and  so  that  the  mag- 
nitudes could  be  related  to  particle  size  through  Mie 
theory  predictions.  In  order  to  relate  the  detected  signal 
amplitudes  to  absolute  scattered  power,  it  was  necessary  to 
calibrate  the  photomultiplier  tube  assembly  with  regard  to 
absolute  anode  sensitivity.  In  order  to  relate  the  c-  . ^i- 
cal  levels  to  photoelectron  rate,  it  was  necessary  to  measure 
the  product  of  quantum  efficiency  and  collection  efficiency 
of  the  PMT . These  measurements  have  already  been  discussed. 
Extreme  care  was  used  in  eliminating  low-level  ground  loops, 
room  light  effects  and  DC  signal  drifts.  The  measurement 
was  made  with  very  low  velocity  (1.0  m/sec)  and  filter  band- 
width chosen  to  pass  only  the  pedastal.  This  was  done  after 
observation  of  the  signals  at  wider  bandwidth  showed  that 
they  had  good  fringe  visibility.  The  extreme  low-pass  fil- 
tering was  accomplished  using  the  filter  in  a TEKTRONIX  1A7A 
plug-in  unit  as  illustrated  in  Figure  10.  The  sweep  speed 
was  made  just  fast  enough  to  nearly  fill  the  display  with  a 
single  burst  so  that  few  were  missed.  (The  bursts  remained 
non- overlapping,  except  in  rare  events,  down  to  the  lowest 
levels  measurable  (5  p volts) . 

★ 

Sec  "Scaling"  in  Appendix  B. 
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Figure  11  summarizes  the  results  of  cumulative  rate/ 
amplitude  measurements  obtained  using  the  a rgon  back-scatter 
system  just  described.  The  procedure  included  using  low 
velocity  and  low  filter  bandwidth  to  obtain  amplitude  meas- 
urements from  particles  which  would  only  scatter  a few  pho- 
tons during  the  short  transit  time  at  typical  transonic 
velocities.  Ten-second  counts  were  made  at  each  trigger 
threshold  settings.  Three  consecutive  runs  made  over  a 
period  of  l*j  hours  showed  excellent  agreement  and  repeat- 
ability. 

Particle  size  computations  were  made  based  on  the 
n * 1.55  Van  DeHulst  curves.  The  parallel  polarization 
values  at  166°  scattering  angle  were  used  to  be  consistent 
with  the  experiment.  Also,  calibrated  PMT  and  optical  sys- 
tem parameters  were  used.  The  right  hand  bars  in  Figure  11 

are  the  calculated  values  for  a particle  through  the  center 

2 

of  the  probe  volume.  The  left  hand  bars  are  1/e  values  of 
amplitude  that  show  the  minimum  peak  amplitude  that  the  same 
particle  size  could  produce.  (The  collecting  optics  pinhole 
was  matched  to  the  l/e^  pulse  diameter.) 

5.1.5  Scaling  of  Burst  Amplitude  Data 

The  optical  parameters  of  the  laboratory  burst  amplitude 
(aj ) distribution  measurement  were  chosen  as  closely  ao  pos- 
sible to  coincide  with  future  4T  transonic  wind  tunnel  optics 
after  discussion  and  with  the  help  of  Mr.  Don  Brayton.  The 
wavelength,  beam  interi  ction  angle,  and  focal  beam  diameters 
were  nominal.  Velocity,  laser  power,  and  collecting  f/no  must 
be  scaled.  The  two  largest  unpredictables  are  the  difference 
between  the  natural  aerosol  in  the  lab  test  and  the  transonic 
test  and  the  differences  in  on-axis  backscatter  and  small- 
angle  off-axis  backscatter  v?ith  regard  to  probe  volume  size 
and  definition.  Neglecting  the  last  two  factors  (although 
they  are  significant)  we  can  scale  velocity,  laser  power,  and 
collecting  f/no  to  obtain  a typical  data  set. 
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Figure  11.  Burst  Amplitude  Data  for  Low-Speed 

Natural  Air  (1.8m/sec)  and  Dita  Scaled 
to  32S  m/sec. 
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The  assumed  scaling  parameters  are  given  in  Table  7. 


f • Lab 

• Velocity  1.8m/sec 

i Signal  Frequency  143  KHz 

Signal  Period  7 sec 

Burst  Duration  175  sec 

Collecting  f/no  3.7 

* Laser  Power  150mw 

Quantum  Efficiency*  0.04 


4T 

325m/sec 
25.8  MHz 
39naec 
0.969  sec 
6 

2w 

0.12 


Ratio 

180.6 

180.6 

5.54  x 10~5 
5.54  x 10'J 
0.62 
13.3 
3 


Table  7.  Scaling  Burst  Amplitude  Data. 


The  procedure  is  as  follows t The  data  rate  at  a given 
value  of  aj  is  scaled  by  the  velocity  (180).  The  value  of 
aj  at  a given  rate  is  scaled  by  the  laser  power  and  the  square 
of  the  f/no  ratio  (5.1).  Thus  the  log-log  plot  is  the  same 

but  translated  up  by  the  velocity  and  translated  to  the  right 

by  the  power  ratio.  The  particle  size  bars  are  translated 

in  the  same  manner.  The  result  of  these  operations  is  illus- 

trated by  an  additional  line  in  Figure  11 . 

The  data  given  in  Figure  11  in  load  voltage  may  be 
converted  to  the  more  significant  quantity  Xj  the  peak 
pedastal  photoelectron  rate  in  photoelectrons/sec  by  the 
formula 

Xj"  hv  (50n)S  (5'1) 

A 

where  for  the  tube  used  (see  PMR-20,  Table  5)  we  obtain 
X j “ 2 . 0 x lO^a  j . Figure  11  indicates  the  locations  of  the 
X » 108/sec  rate  (maximum  for  state-of-the-art  photon-count- 
ing) and  also  a .line  corresponding  to  the  lower  limits  of 
burst  counters  with  a 25  MHz  signal  frequency. 


including  collecting  efficiency.  This  assumes  use  of 
RCA  8850  in  4T. 
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5.2  BACKGROUND  AND  SIGNAL  LEVEL  EXPERIMENTS 

In  April  of  1975  a trip  was  made  to  the  Arnold  Center 
to  check  out  the  Ortec  and  Systron  Donner  equipment  and  the 
photon  counter  portion  of  the  experimental  processor.  In 
addition  we  wished  to  observe  photon  signals  from  a back- 
scatter  LV  system  which  Virgil  Cline  and  others  at  the 
Arnold  Center  had  constructed  for  use  in  the  IT  transonic 
wind  tunnel.  The  maximum  photon  discriminator  pulse  rate 
is  just  over  100  MHz  so  it  was  important  to  see  what  the 
rates  were  from  a typical  real  LV  system. 

5.2.1  Equipment  Checkout 

The  Wavetek  oscillator  (WOS)  was  used  to  drive  the 
Systron  Donner  Counter  (SDC) . The  SDC  operated  in  linear 
agreement  with  the  dial  reading  of  the  WOS  to  230  MHz.  The 
WOS  was  applied  to  the  amplifier  discriminator  with  the  SDC 
counting  the  discriminator  output  pulBes.  This  produced 
reliable  behavior  to  100  MHz  with  low  level  input  (unrecorded); 
with  increased  amplitude,  the  range  was  extended  to  125  MHz. 

The  RCA  8644  and  8850  tubes  were  not  used  during  this 
trip.  Virgil  Cline's  set-up  included  an  RCA  931A  so  this 
tube  was  used  exclusively.  A variety  of  arrangements  of  50  U 
cables,  attenuators,  terminators,  cable  lengths,  and  ampli- 
fier combinations  were  tried  out  to  minimize  after-pulse 
ringing  which  was  at  first  a serious  problem.  The  bandwidth 
restriction  and  pulse  shaping  effects  of  the  9302  amplifier 
increase  the  pulse  width  slightly  from  that  amplified  and 
delivered  by  the  9301  preamp,  but  also  reduces  the  high- 
frequency  ringing. 

The  effect  of  varying  the  discriminator  threshold  was 
observed.  This  experiment  dramatically  demonstrated  the 
significance  of  the  pulse  height  variation.  When  the 
threshold  was  set  low  (50-100  mv  after  lOx  preamp  and  20x 
amp) , the  oscilloscope  revealed  secondary  firings  of  the 
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discriminator  at  the  spacing  of  the  first  round-trip  pulse 
reflection  - 20  nsec  (even  with  all  the  precautions  of  pre- 
cision cable  and  terminations) . The  secondary  pulses  were 
caused  by  the  large  single-photo-electron  pulses  whose 
reflections  were  also  above  the  threshold.  These  were 
removed  successfully  by  setting  the  threshold  up  to  400 
mv.  The  counted  photon  rate  (disc  out  to  SDC)  decreased 
by  a factor  of  20. 

The  result  of  the  threshold  experiment  shows  that 

(a)  The  931A  cannot  be  efficiencly  used  for  fast 
photon-counting  work  because  pulse-height 
variation  forces  one  to  set  the  threshold  so 
high  that  most  of  the  signal  pulseB  are  missed. 

(b)  The  experiment  of  comparing  the  8650  with  other 
tubes  for  count  rate  cannot  be  performed  with 
any  meaning  until  criteria  for  secondary 
threshold  crossings  are  defined.  In  photon  cor- 
relation work  this  criteria  is  expected  to  be 

so  stringent  that  the  effective  useable  quantum 
efficiency  of  the  8850  will  be  many  times  greater 
than  the  931A. 

ORTEC  9 specifies  nsec  as  the  typical  dead-time  for  the  9302 
amp/discriminator  and  guarantees  10  nsec  maximum.  Driving 
the  input  with  an  oscillator  and  getting  output  pulses  at 
125  MHz  would  seem  to  indicate  an  8 nsec  dead  time  for  the 
unit  we  obtained.  However,  when  the  discriminator  output 
was  observed  on  the  oscilloscope  at  hiqh  light  levels 
(X  ■ 50  x 10 6 ) the  minimum  separation  observed  between  the 
triggering  pulse  and  the  next  one  was  9 to  10  nsec.  We 
c nsider  the  combined  experiments  to  verify  the  "typical” 
number  of  9 nsec. 

Check  out  of  the  photon  counter  indicated  that  it  was 
only  performing  satisfactorily  to  count  rates  of  about  60 
MHz.  The  circuitry  was  later  redesigned  and  reconstructed 
to  operate  to  over  100  MHz  count  rates. 
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5.2.2  Background  Rates 

This  was  a significant  experiment.  We  wished  to 
determine  background  rates  that  could  be  expected  with  a 
practical  LV  optical  set-up  in  use  at  the  A£DC.  With  50 
mv  discriminator  threshold  (RCA  931A  with  200x  amplification 
before  the  threshold)  we  obtained  the  following  results  by 
counting  the  discriminator  pulses  with  the  Systron  Donner 
counter . 

Condition  Rate 

1.  PMT  covered  (Dark)  0.10  MHz 


2.  Room  lights  out  and  black 
curtain  in  front  of  window 

3.  Room  lights  off,  window  light 
on  wall  in  field  of  view 

4.  Same  ac  (3)  but  with  laser  on 
(signal  from  natural  air) 

5.  Same  as  (3)  but  with  over- 
head room  lights  on 


0.15  MHz 


1.0  to  2.5  MHz 


« 25  MHz 


“ 50  MHz 


6.  Same  as  (4)  but  with  typical 
plexiglass  window  inserted  in 
beams.  (No  noticeable  change.)  “ 25  MHz 

The  signalB  obtained  are  in  an  acceptable  range  of  rates  for 
photon  counting.  The  exception  was  the  background  count 
with  the  room  lights  on.  This  experiment  shows  the  need 
for  careful  shielding  of  the  PMT  optics. 

5.3  FEASIBILITY  DEMONSTRATION  OF  THE  DUAL  CORRELATE  AND 
SUBTRACT  MODE 

On  September  2,  1975  a trip  was  made  to  the  Arnold 
Center  with  the  experimental  photon  processing  system, 
pj  Although  deterministic  laboratory  checkout  had  indicated 

that  all  the  developmental  bugs  had  been  eliminated  from 
the  circuitry,  this  was  not  the  case.  Additional  debugging 
of  the  circuits  was  required.  No  measurements  were  obtained 
until  the  last  half  of  the  day  September  5,  1975.  The  system 
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was  used  that  day  to  make  mean  flow  measurements  In  a 
small  laboratory  jet  under  weak  signal  conditions  where 
the  ARO  Model  8 burst  counter  processor  would  not  make  any 
measurements.  A small  amount  of  data  was  obtained  for  bias 
error  and  variability  error  evaluation.  The  mean  velocity 
estimate  differed  in  one  instance  by  4%  from  the  velocity 
measured  by  the  ARO  counter  at  3 times  greater  laser  power. 
This  is  good  agreement  since  the  ARO  processor  was  barely 
producing  data  (1  data  point  validated  per  3 seconds  at  the 
higher  power) . 

The  variability  error  agreed  very  well  with  the  vari- 
ability theory  which  is  discussed  in  section  6.0.  The  experi- 
ment is  described  in  more  detail  below. 

5.3.1  The  Experiment 

An  off-axis  backecatter  (fringe-type)  LV  system  opera- 
ting with  a wavelength  of  514.5  nm  and  fringe  spacing  of 
6.9  micrometers  was  assembled  for  a laboratory  test  with  a 
pipe  jet  air  flow  (20  ft/sec,  mean  velocity).  The  shop  air 
jet  appeared  visually  to  scatter  less  light  than  the  room 
air  so  that  the  beam  crossing  location  was  nearly  invisible 
in  comparison  with  the  other  beam  locations.  After  align- 
ment, unusually  weak  signals  were  obtained  from  an  RCA  931A 
PMT  even  with  1.7  watts  of  laser  power.  When  a low-pass 
filter  (1  MHz)  was  used  to  increases  the  reliability  of  trig- 
gering on  signal  bursts  and  when  the  room  lights  were  extin- 
guished,* the  validated  data  rate  from  the  ARO  model  8 
processor  (3%  window)  was  approximately  1 per  3 sec.  The 
Doppler  frequency  was  slightly  less  than  3.  MHz.  The  level 


The  signal  levels  were  sufficiently  low  that  room  light 
was  the  predominant  source  and  neither  type  of  processor 
worked  properly  with  them  on  (mean  photon  count  " 100  times 
higher  with  room  lights) . 
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of  turbulence  was  approximately  7»8%  as  measured  by  the 
ARO  Model  8 at  * 1.7  watts. 

The  demonstration  of  concept  feasibility  consisted  of 
obtaining  mean  velocity  measurements  with  the  laser  power 
turned  down  to  540  mw  and  100  mw  (corresponding  to  mean 
photo  electron  rates  of  547  KHz  and  105  KHz.)  Under  the 
former  condition,  1 validated  Model  8 measurement  was 
obtained  in  30  sec.  and  at  100  mw,  none  were  obtained  in  2 
minutes  of  waiting  and  twiddling  with  the  trigger  threshold. 
The  photon  processor  was  able  to  make  a mean  velocity  esti- 
mate at  both  100  mw  and  500  mw  by  operating  it  in  a feed- 
back control  loop  mode  with  a human  operator  adjusting  the 
Wavetek  oscillator  over  a period  of  a few  seconds  while 
watching  the  accumulator  display  lights. 

5.3.2  Statistical  Error 


A bias  error  occurs  when  the  long  term  average  value  of 
a measurement  differs  from  the  true  average  of  the  quantity 
being  measured.  The  signal  levels  were  so  low  that  insuf- 
ficient data  was  collected  to  establish  very  accurate  long 
time  averages  for  either  the  Model  8 processor  data  or  the 
photon  processor.  Also,  there  is  no  reason  that  the  ARO 
processor  was  exactly  correct  (3%  error  window,  and  very 
low  signal  levels).  In  spite  of  the  above,  we  have  compared 
the  30  sec  (10  measurement)  average  velocity  obtained  from 
the  ARO  Model  8 at  1.7  wattB  of  laser  power  with  the  values 
obtained  by  human  feedback  adjustment  of  the  variable  oscilla- 
tor. The  values  agree  within  4%  for  the  540  mw  case  and 
within  20%  for  the  100  mw  case. 

A much  better  statistical  measure  of  the  variability 
error  was  obtained  for  comparison  with  theory.  With  the 
oscillator  frequency  set  at  the  value  which  nulled  the 
accumulator,  a scries  of  10  one-second  measurements  were 
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made  for  each  laser  power  level.  The  value  accumulated  each 
time  was  recorded/  and  the  sample  rms  deviation  was  com- 
puted as  shown  in  Table  7. 

Section  6.0  provides  a theoretical  equation  for  the 
rms  deviation  of  the  accumulator  results  as  follows t 


/5n  <x>  At 


(5.2) 


N ■ T/Ar is  the  total  number  of  clock  periods  during  the 
measurement  time.  This  simplified  theory  assumos  a con- 
stant amplitude  light  source  and  is  therefore  only  approxi* 
mate.  We  have  obtained  the  following  comparisons: 


rl 

100  mw 


500  mw 


o Theory 
50.2 
285 


a Experiment 
48.4 
309.0 


Table  7.  Comparison  of  Theoretical  and  Experimental 
Standard  Deviation  of  Accumulator  Values. 


5.4  COUNTER  BIAS  AND  PHOTON  CORRELATION 


When  the  feasibility  demonstration  experiment  just 
described  was  performed,  there  were  several  minor  opera- 
tional problems  such  as  overheating  and  occasional  inter- 
mittent accumulator  operation  which  occurred  at  clock  fre- 
quencies higher  than  those  required  for  the  demonstration. 
These  problems  were  fixed  later  by  minor  circuit  changes  and 
the  addition  of  more  powerful  cooling  fans.  On  October  10 
and  13,  1975  SAI  and  GW  Electronics  personnel  performed  addi- 
tional laboratory  measurements  to  verify  proper  operation  of 
all  system  functions  at  speeds  up  to  the  guaranteed  minimum 
clock  frequency  of  40  MHz.  Two  significant  additional 
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Clock  Frequency  ~ MHz 

Figure  12.  Efficiency  of  Four-Bit  Counter  Vs.  Clock  Frequency  and  Data  Rate 
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measurements  were  conducted  using  an  RCA  8645  PMT  with 
the  ORTEC  9301  preamp,  and  a light  emitting  diode  source 
(LED) . 

5.4.1  Counter  Bias  Vs.  Clock  Frequency 

The  design  of  the  4-bit  counter  circuit  uses  two 
alternating  counters  with  a small  dead  time  between  count 
intervals.  In  order  to  measure  this  effect  we  placed  the 
processor  in  MEAN  mode,  so  that  it  merely  accumulated  the 
sum  of  all  the  n^  during  the  counting  interval,  and  com- 
pared the  sum  with  the  total  accumulated  by  the  Systron  Donner 
laboratory  counter  (SDC)  during  the  same  interval.  We  did 
this  for  several  different  mean  photoelectron  rates  and  sev- 
eral different  system  clock  frequencies  and  parametrically 
plotted  the  ratio  of  the  processor  total  to  the  SDC  total. 

This  data  is  illustrated  in  Figure  12.  The  figure  reveals 
three  significant  effects  upon  close  scrutiny.  First,  over 
most  of  the  data,  the  count  efficiency  i3  not  a function  of 
photon  rate.  Second,  at  the  very  low  clock  frequencies  and 
high  photon  rates,  the  count  saturates  at  15  (4  bit  counter) 
and  data  is  lost.  Third  the  slight  high  frequency  slump  is 
predicted  almost  exactly  by  assuming  an  interperiod  dead  time 
of  2.3  nsec. 

5.4.2  Photon  Correlation 

The  LED  was  modulated  at  3.000  MHz  with  the  bias  dc 
level  adjusted  to  give  a ratio  of  the  peak-to-peak  ac  light 
modulation  to  steady  light  ratio*  of  0.18.  The  LED  was 
moved  away  from  the  shielded  PMT  (with  pinhole)  until  the 
mean  photo-election  count  rate  was  5x10^  or  1.67  photoelec- 
trons per  cycle.  The  photon  processor  was  sequentially 


As  measured  by  an  oscilloscope  observation  of  a rela- 
tively high  level  signal  output  from  the  PMT  directly  to  a 
load  resistor. 
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programmed  through  all  delay  values  with  the  processor  in 
AUTOCORRELATE  MODE,  the  system  clock  at  40,00  MHz,  and  the 
measurement  interval  per  delay  of  0.1  sec.  The  total  photo- 
electron count  was  measured  by  the  SDC  simultaneously  for 
each  delay  value.  Since  the  mean  photoelectron  rate  drifted 
during  the  time  the  dat-a  were  collected,  the  autocorrelation 
results  were  normalized  by  dividing  by  the  square  of  the 
totalized  photoeiectron  count  for  each  delay  interval. 

The  results  of  this  experiment  are  plotted  in  Figure  13 . 
It  is  significant  that  with  an  average  of  1.7  photoeiectron ic/ 
cycle,  the  entire  measurement  would  have  required  only  2.0 
seconds  if  it  had  been  sequenced  by  a computer.  The  value 
plotted  at  500  nsec  delay  is  actually  3 separate  0.1  sec 
runs.  This  is  an  indication  of  the  low  variability  error 
for  this  particular  experiment. 
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6.0  DISCUSSION 

In  this  section  we  discuss  photomultiplier  tube 
effects,  burst  counter  threshold,  and  errors  in  the  meas 
urement  of  mean  flow  with  the  new  Dual  Correlate  end  Sub 
tract  technique. 


6.1  PHOTOMULTIPLIER  TUBE  NOISE  EFFECTS 

There  have  been  so  many  written  statements  to  the 
effect  that  photomultiplier  tubes  are  "nearly  noiseless" 
that  we  may  begin  to  believe  it  doesn't  matter  much  which 
tube  is  used  except  for  speed,  size,  ruggedness,  quantum 
efficiency,  and  cost.  For  all  low  level  LV  signal  detection 
schemes  for  photon-limited  and  photon-resolved  signals 
nothing  could  be  less  true.  The  significant  factors  in  addi- 
tion to  the  above  mentioned  are  dynode  collection  efficiency, 
dynode  pulse  charge  gain  statistics,  pulse  transients  (ring- 
ing) , correlated  after  pulsing,  and  total  dynode  gain.  The 
effects  of  pulse  gain  statistics  and  dynode  collection  effi- 
ciency may  be  serious  even  for  signals  in  the  Gaussian 
regime  since  they  may  easily  reduce  the  signal  power  to  noise 
ratio  (see  equation  2.13  for  definition)  by  a factor  of  4 or 
more.  References  [16-183  should  be  read  for  thorough  dis- 
cussion of  effects  to  be  considered.  In  light  of  the  fol- 
lowing, we  recommend  that  the  newer  tubes  with  small  pulse 
height  variation  and  high  collection  efficiency  be  used  in 
future  LV  systems. 

6.1.1  Dynode  Collection  Efficiency 

The  effective  quantum  efficiency  n for  photon  counting 
is  the  product  of  nc,  the  cathode  efficiency  usually  considered 
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and  F,  the  dynode  pulse  efficiency.  This  dynode  efficiency 
F is  less  than  unity  because  some  of  the  cathode  photo- 
electrons are  collected  by  the  dynodes  without  producing  an 
output  pulse.  The  point  missed  by  many  investigators  is 
that  the  product  ncF  is  the  correct  quantity  to  use  in  SNR 
equations  (or  in  the  computation  of  X(t))  even  at  high 
signal  levels.  It  is  not  difficult  to  measure  nc  by  direct 
cathode  sensitivity  measurements.  The  quantity  F is  not 
specified  in  typical  manufacturers  specifications.  Our 
experience  (see  Section  5.1.1)  is  that  n may  be  less  than 
half  of  the  specified  quantum  efficiency  due  to  either 
cathode  fatigue  and/or  poor  dynode  collection  efficiency. 

6.1.2  Single  Photoelectron  Gain  Variation 

The  charge  delivered  to  the  anode  by  the  ith  photo- 
electron pulse  is  e g^,  where  g^  is  the  random  charge  gain. 
The  relationship  between  the  DC  current  gain  G for  a given 
tube  and  the  average  value  <g^>  of  the  single  photoelectron 
pulse  gain  may  be  shown  to  be 

<g^>F  =*  G (6.1) 

The  Bhape  of  probability  density  for  determines 
several  significant  things  including  the  value  of  the  pre- 
multiplier term  in  the  signal-to-noise  ratio  equation  2.13. 
Reference  [14’J  pages  65  and  66  wou-d  imply  that  for  a typical 

copper  berylium  dynode  tube  with  stage  gains  of  3-4  the 

2 2 

reduction  factor  <g>  /<g  > will  be  in  the  range  0.5  to  0,76. 
This  is  not  an  insignificant  factor  and  should  be  included 
carefully  in  PMT  selection.  For  some  of  the  newer  tubes 
with  high  first-stage  gain,  the  factor  would  be  nearly 
unity . 
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For  very  low-level  signal  detection  and  photon  count- 
ing, the  pulse  height  variation  is  harmful  in  several  ways. 
First,  the  variation  of  height  can  cause  false  triggering 
on  the  anomalously  large  pulses  in  a low-level  burst- 
counter  or  transient  recorder  system.  In  a photon  counting 
system  a timing  error  occurs  because  of  threshold  crossing 
at  different  portions  of  the  pulse  waveform.  Another  consid- 
eration is  the  fact  that  for  a low-gain  PMT  the  Gaussian 
electronic  preamplifier  noise  may  easily  be  the  same  order 
of  magnitude  as  the  smaller  single  photoelectron  pulses. 

This  is  true,  for  example,  with  the  RCA  8644  or  8645  tubes. 
Another  problem  with  large  pulse  height  variation  is  that  in 
very  fast  system,  the  pulses  have  a ringing  tail  no  matter 
how  careful  we  are  in  the  physical  arrangement.  When  we  set 
a discriminator  threshold  above  the  ringing  tail  of  the 
larger  pulses,  we  then  miss  the  smaller  oneB.  Photon  count- 
ing systems  thus  need  great  care  in  circuit  design  and 
construction  details  and  a tube  with  Bmall  pulse  height 
variation. 

6.2  BURST  HEIGHT  DISTRIBUTION  EFFECTS 

It  is  clear  from  the  experimental  data  shown  in  Figure 
11  on  burst  amplitude  vs.  rate  of  occurrence  that  with 
natural  aerosols  there  may  be  a few  large  signals  with  ade- 
quate "signai -to-noise-ratio"  for  burst-counter  operation 
and  many  smaller  bursts  which  produce  only  a few  photo- 
electron pulses.  There  are  the  questions  of  how  do  we  com- 
pare a photon  counter  processor  and  a burst  counter  processor, 
how  do  we  determine  the  source  distribution  of  the  signals 
which  are  rejected  by  each  type  of  processor  and  which  are 
accepted . 
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Appendix  B addresses  itself  to  the  above  questions 
and  produces  some  answers  and  some  more  questions.  For 
example,  neglecting  PMT  pulse  pile  up  and  discriminator 
saturation  and  dead  time  effects,  the  appendix  shows  that 
for  the  experimental  data,  the  slope  of  the  amplitude  rate 
distribution  is  sufficiently  steep  that  the  accumulator 
sum  places  emphasis  on  the  smaller  scatterers,  which  produce 
fewer  photons.  The  appendix  does  not  answer  the  question 
of  what  happens  to  a photon  counting  processor  when  the 
burst  iB  photon  limited  and  not  large  enough  to  go  above 
the  discriminator  threshold  and  stay  (in  which  case  it 
would  contribute  negligibly  to  the  result.)  We  expect  a 
phase  reversal  contribution  in  some  cases  where  the  valleys 
of  the  signal  produce  threshold  crossings  by  single  photo- 
electron  pulses  while  the  peaks  remain  above  threshold  and 
produce  no  pulses.  For  cases  where  the  rate  amplitude 
distribution  has  sufficiently  negative  slope  in  the  photon- 
limited  regimes  and  below,  we  would  expect  the  as  yet  unpre- 
dictable effects  to  be  negligible  in  the  sum,  but  obviously 
further  work  is  needed. 

In  contrast  to  the  maximum  photoelectron  rate  which 
a photon  counting  system  can  utilize  due  to  electronic  limi- 
tations, the  required  rate  appears  to  be  much  less.  The 
scaled  data  in  Figure  11  indicates  a data  rate  of  approxi- 
mately 100,000  signal  burst/sec  at  the  peak  rate  level  of 
_>  107  photo-electrons/second . Since  the  scaled  burst  dura- 
tion is  1 ysecond,  this  corresponds  to  a duty  factor  of  0.1 
at  this  level.  These  figures  are  similar  to  the  parameter 
values  chosen  in  our  "typical"  variability  error  example  in 
section  6.3  where  the  preliminary  theory  indicates  1%  rms 
variability  error  in  mean  flow  measurement  with  0.5  seconds 
of  data.  Our  present  analysis  indicates  that  with  the  exam- 
ple optical  system  the  contributing  particle  size  would  be 
primarily  in  the  0.2  - 0.3u  particle  size  range. 
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The  lower  threshold  of  burst  counter  processors  is 

O 

shown  as  approximately  4 mv  or  8 x 10  photoelectron  pulses 
sec  peak  pedestal  in  Figure  11.  This  threshold  was  origi- 
nally arbitrarily  based  on  the  assumption  of  needing  4 

2 

photo  electrons/cycle  at  the  1/e  pedastal  locations  with 

a 25  MHz  signal)  i.e.  we  assumed  the  need  for  greater  than 

4 photoelectron  pulses/cycle  over  the  1/e  burst  width. 

2 8 8 

This  leads  to  a peak  pedastal  value  of  X^-e  x 10  « 7.4x10  . 

Independently  of  the  above  assumptions , Franz  Durst 

and  others  in  discussions  at  the  1975  Minnesota  Conference 

[19]  indicated  that  the  lower  limit  of  reliable  burst-counter 

operation  was  at  average  signal-to-noise  power  (during  peak 

signal)  ratios  of  10.  From  equation  2.13  with  0.707  assumed 
2 2 

for  <g^>  /<g^>,  we  would  obtain  a requirement  for 
X,  - 10  x /I  ■ 14.1  as  a threshold.  This  would  indicate 
approximately  7-8  mv  pedastal  signals  for  the  RCA  8644  tube 
used  would  be  required  and  our  scaled  estimate  of  300  accept- 
able burst/sec  as  too  high. 

In  reality,  there  is  no  absolute  "threshold*'  for 
acceptable  burst  counter  operation.  The  probability  of  cor- 
rect detection  is  a continuous  function  of  signal  amplitude, 
background  level,  and  low-level  signal  distribution.  The 
problem  of  analysis  for  burst  counters  is  formulated  in  more 
detail  in  Appendix  B,  but  no  attempt  at  its  solution  is  pre- 
sented. Further  work  in  this  area  both  analytically  with 
simulation  is  recommended. 

6 . 3 STATISTICAL  ERRORS 

The  two  principle  types  of  error  which  arise  in  statis- 
tical measurements  are  bias  error  and  variability  error. 

Bias  error  is  a term  which  refers  to  the  difference  between 
the  statistical  expectation  of  the  measurement  system  output 
and  the  desired  average  value  being  measured.  The  variability 
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error  is  the  rms  value  of  the  random  deviation  of  a 
specific  experimental  result  from  the  statistically  expected 
value.  For  ergodic  random  processes,  the  variability 
error  converges  to  zero  in  the  limit  of  infinitely  long  data 
collection  time;  but  it  converges  to  an  acceptable  level 
(which  must  be  defined  in  the  measurement  objectives)  within 
a finite  measurement  time.  The  bias  error  cannot  be  removed 
by  further  averaging  but  it  can  often  be  removed  by  analyti- 
cal compensation  or  by  experimental  calibration  when  it  is 
3mall  compared  with  the  desired  quantity.  In  general, 
analysis  of  both  types  of  errors  is  required  in  any  statis- 
tical measurement.  Bias  errors  are  evaluated  by  using  more 
precise  signal  and  system  models  in  determining  expected 
values  of  the  output. 

Error-analysis  of  a Dual  Correlate  and  Subtract 
processor  will  be  presented  in  more  detail  in  the  final 
report  for  Contract  NASAl-13737,  between  NASA  Langley 
Research  Center  and  Science  Applications,  Inc.,  which  is  not 
yet  complete.  In  this  section  we  briefly  indicate  the  nature 
of  bias  errors  which  must  be  investigated  and  summarize  pre- 
liminary variability  error  analysis  results  which  have  been 
obtained  for  NASA. 

6.3.1  Bias  Errors 

There  are  several  potential  sources  of  bias  error  for 
the  processor  we  have  designed.  First,  the  value  of  n^  does 
not  represent  the  number  of  photoelectrons  detected  at  time 
kAx.  The  detection  events  occur  in  the  time  interval 
(k-1) Ax<t<kAx . The  effect  of  this  time  discretization  is 
similar  to  that  of  our  previous  work  [20D  in  that  the  expected 
value  of  an  autocorrelation  estimate  at  a discrete  delay  is 
a weighted  average  of  the  true  autocorrelation  function  over 
an  interval  about  that  point.  The  weighting  function  is 
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triangular,  extending  ±Ax,  as  is  shown  in  equation  2.24. 

The  effect  of  a triangular  weighting  function  on  the  cor- 
relation  estimate  is  that  of  a [sin  (irAif )/irArf]  low-pass 
filter.  For  Ax  ■ l/4fm  the  cosinusoidal  portion  of  the 
autocorrelation  estimate  is  thus  attenuated  by  a factor  of 
0.81  and  lower  signal  frequencies  are  attenuated  less  than 
higher  ones.  A small  bias  of  the  mean  frequency  estimate 
will  thus  be  developed  which  vanishes  for  small  turbulence 
intensities  or  for  smaller  values  of  Ax . 

Additional  bias  error  will  be  induced  by  time  jitter 
in  the  detection  of  photoelectron  pulse3.  For  a fixed 
threshold  trigger  circuit,  different  portions  of  the  pulse 
leading  edge  will  trigger  the  discriminator  if  the  amplitude 
of  the  single  photoelectron  pulses  varies.  With  the  presently 
available  PMT's,  this  jitter  can  be  no  more  than  the  rise 
time  of  less  than  2 nsec  but  can  be  reduced  to  a frac- 
tion of  a nsec  with  newer  tubes  with  small  pulse  height 
variation.  The  effect  is  small  compared  with  the  Ax  effect 
even  at  the  fastest  range  of  10  nsec. 

The  threshold  discriminator  dead  time  will  produce 
bias  errors.  It  is  obvious  that  if  the  dead  time  were  20 
nsec,  there  would  be  no  lag  products  at  a 10  nBec  delay 
interval.  The  functional  dependence  of  the  bias  error  has 
not  been  derived  for  this  system  yet.  On  the  highest  range, 
with  Ax  - 10  nsec,  signal  frequency  ■ 25  MHz  and  with  a dead 
time  of  9 nsec,  the  estimate  will  be  distorted  because 
delay  1 implies  events  in  the  range  (0.0  Ar , 2.0  Ax)  while 
delay  3 refers  to  events  separated  in  time  by  (2.0  Ax,  4.0Ax). 
Clearly,  the  delay  equal  1 term  will  be  reduced  appreciably 
on  the  faster  ranges.  This  effect  may  be  seen  in  Figure  13 
where  the  first  delay  term  appears  slightly  low.  It  is 
possible  that  correction  curves  may  be  derived  in  conjunc- 
tion with  experimental  calibration  results,  but  further 
work  is  needud  to  assess  the  limitations. 
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"In  the  hardware  design  we  have  included  a provision 
for  changing  the  relative  width  of  At  with  respect  to  the 
signal  frequency.  This  may  be  accomplished  simply  by 
increasing  the  numbers  of  clock  delays  for  the  indices  p 
and  q while  keeping  their  ratios  constant.  For  example, 
with  At  ■ 1/4  period  in  the  dual  correlate  mode,  p~3,  q«l> 
by  changing  p«6,  q ■ 2 At  becomes  1/8  period  at  the  mean 
signal  frequency.  The  price  for  the  change  is  the  reduc- 
tion of  the  maximum  potential  signal  frequency  from  25  MHz 
to  12.5  MHz  and  increased  variability  error,  but  bias  error 
effects  due  to  At  width,  after-pulse  correlation,  and  dis- 
criminator dead  time  may  all  be  significantly  reduced. 

The  analytical  formula  for  the  statistical  expecta- 
tion of  the  output  of  s photon  correlator  was  obtained 
using  3.7: 


<COS  (0  j T > 


p^(u)  COS  U)T  do) 


cos  iii  t (6.3) 
m 


The  specific  result  depended  on  both  the  assumption  that 
Pu(u)  was  Gaussian  and  the  assumption  that  particle  rate 
bias  effects  did  not  keep  from  being  the  same  shape. 

It  is  not  difficult  to  show  that  cos  u>  t factor  will  always 
result  if  p(u>)  is  symmetrical  about  its  mean  tu  . The  ques- 
tion of  concern  is  "what  happens  if  p(ui)  is  skewed?"  Pre- 
liminary efforts  indicate  that  there  are  error  contributions 
from  higher  order  odd  central  momenta  of  p(iu),  whose  magni- 
cude  depends  on  the  deviation  of  the  discriminator  charac- 
teristic from  linearity  and  which  vanish  with  decreasing 
turbulence  intensity  as  expected.  This  work  has  not  been 
refined  to  the  point  of  formal  reporting.  We  feel  intui- 
tively that  the  results  will  show  that  the  problem  is  minor 
for  turbulence  levels  on  the  order  of  10%  or  loss.  For 
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higher  turbulence  levels,  up  shifting  with  a Bragg  cell  or 
other  means  may  be  required  to  keep  p (<ju ) confined  to  the 
linear  region  of  the  discriminator  function.  It  may  also 
be  possible  in  further  work  to  identify  configurations  for 
which  the  discriminator  characteristic  is  more  linear. 

6.3.2  Variability  Error 

In  an  earlier  very  brief  contract  with  NASA  Langley 
(NAS1-13140)  a simplified  variability  error  analysis  was 
performed  on  the  accumulator  sum  under  the  assumption  that 
the  classical  optical  power  was  steady  at  its  mean  value. 

In  that  case  analysis  of  the  variance  led  to  the  simple  rms 
result i 


■ /2H  < X > Ax  (6.4) 

mpq 

which  although  much  oversimplified  has  shown  good  agreement 
with  experiment  (see  Section  5.3.2) . In  the  reBt  of  this 
subsection  we  paraphrase  the  NASA1-13140  report*  in  which 
we  applied  equation  6.4  to  obtain  an  estimate  of  the  frac- 
tional rms  velocity  error,  ■ a^/  U.  The  resulting  approxi- 
mation waB 


E 


U 


+ 1] 


7T  X /TAl 

Ol 


(6.5) 


where 


X 

X 

X 


b 

s 

a 


T 


« background  photoelectron  rate  (6.6) 

* mean  signal  rate  « R<X^>/jT  a 

= < X j > * mean  peak  pedastal  photoelectron  rate 

* NAr  ■ data  collection  time 
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In  order  to  provide  further  intuitive  insight  concern- 
ing equation  6.5  we  include  the  following  example  numbers 
which  might  be  typical  for  a transonic  measurement: 


Xb  - 107/sec  (6.7) 

Xg  - 106/sec 

X - 107/sec 

A 

At  ■ 10~®  sec 

eu  “ °‘01 
T - 0.49  sec 


In  this  example,  the  mean  signal  photoelectron  rate  is  ten 
times  less  than  the  mean  background  photo-electron  rate  and 
is  equal  to  the  average  peak  envelope  rate.  With  this  much 
background  light,  the  assumption  of  constant  P(t)  is  valid 
with  respect  to  the  photon-fluctuation  induced  variability. 

7 

The  selection  of  mean  peak  rate  at  10  means  that  occurrences 

8 

of  photoelectron  count  rates  greater  than  10  /sec  (the  limit 
of  current  hardware  state  of  the  art)  will  be  rare  and  the 
effects  of  nonlinearity  negligible) . The  selected  ratio  of 
X-/X.  ■ 10  implies  that  the  measurement  volume  is  only 
assumed  to  contain  a scatterer  1/10  of  the  time  on  the 
average . 

Even  though  dramatic  improvements  over  the  previous 
example  situation  may  result  from  reduced  background  light, 
the  result  still  indicates  that  practical  measurements  may 
be  obtained  with  a peak  signal  photoelectron  rate  of  107/soc. 
In  order  to  compare  this  with  the  performance  of  a burst- 
counter  system,  we  must  assume  values  for  p,  q and  fm,  tho 
Doppler  frequency.  With  p ■ 3,  and  q ■ 1 the  Doppler  fre- 
quency f is  l/4Ax  ■ 25  MHz.  This  would  result  from  U - 1000 
ft/sec  with  an  optical  sensitivity  of  25  KHz/ft/sec.  This 
peak  photoelectron  rate  assumed  is  thus  0.4  photoelectrons/ 
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cycle  In  the  presence  of  0.4  background  photoelectrons/cycle. 
For  comparison  we  note  that  the  signal-to-noise  ratio  formula 


of  equation  2.14,  given  nP(t)/hv 


(1 


+ cos  u t)  + X.  , 
m d 


would  result  in  peak  average  SNR  (using  2B  = 1/At)  of 

SNR  • 


V 


(6.8) 


In  our  example,  B - 25  MHz  and  X^  - 10  ; if  we  choose  Xfi » 10  , 
a factor  of  100  greater  than  in  our  example,  then  the  SNR  is 
10  at  the  peak  of  the  signal  burst  and  1.35  at  the  1/e  signal 
envelope  points.  Since  this  example  represents  marginal  or 
inadequate  SNR  for  burst  counter  operation  we  deduce  that, 
even  with  100  times  more  scattered  power,  only  the  larger- 
than-average  scatterers  would  contribute. 

Under  conditions  of  less  background  light,  the  burst- 
counter  analysis  would  not  be  improved;  however,  the  photon 
counting  system  results  are  expected  to  improve  considera- 
bly. Thus  we  conclude  that  mean-flow  measurements  with  from 
100  - 1000  times  less  optical  power  are  feasible  with  the 
photon  counting  system. 
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7.0  CONCLUSIONS 

The  effort  under  this  contract  has  resulted  in  several 
significant  conclusions  both  theoretical  and  experimental. 

These  are  summarized  below. 

Theoretical  Poisson  models  have  been  formulated  for 
dual-scatter  LV  signals  with  sufficient  generality  to 
include  all  signal  regimes  from  classical  signals  with 
Gaussian  noise  down  to  photon-resolved  signals  which  require 
photon  counting  techniques. 

Analysis  of  conditional  signal  models  has  provided 
instantaneous  signal-to-noise  ratio  and  other  statistical 
parameter  equations  not  previously  available  in  the  LDV  lit- 
erature . 

Analysis  of  Unconditional  models  has  indicated  the  prac- 
tical feasibility  of  improving  detection  sensitivity  over 
burst-counter  techniques  by  a factor  between  100  and  1000  by 
using  photon  counting  techniques. 

A detailed  analysis,  including  explicit  mention  of 
approximations  required,  of  the  expected  value  of  direct  pho- 
ton correlation  of  LV  signals  from  turbulent  flow  is  provided. 
This  both  supports  and  points  out  limitations  of  previous 
results  given  by  Pike  and  others. 

A new  Dual-Correlate-and-Subtract  technique  for  new  fre- 
quency measurement  which  utilizes  only  two  delay  values  and 
one  real-time  high  speed  multiplier  without  clipping  has  been 
theoretically  shown  to  be  feasible.  This  feasibility  is 
based  on  evaluation  of  the  expectation  for  system  model  to 
show  convergence  to  mean  frequency  estimate,  and  on  rms  vari- 
ability analysis  which  indicates  practical  measurement  times 
on  the  order  of  0.5  seconds  in  transonic  flow. 

An  electronically  switchable,  multipurpose,  high-speed 
electronic  signal  processor  with  the  potential  speed  for  10 
nsec  time  .esolution  and  10  nsec  real  time  4 bit  multiplication 
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has  been  designed  and  constructed  using  slow  ECL  integrated 
circuits  on  special-purpose  100  MHz  wirewrap  panel.  The  sys- 
tem implements  the  dual- correlate-and-subtract  technique. 

It  also  has  the  potential  of  standard  photon  correlation 
without  clipping  by  computer  controlled  scanning  of  the 
delays.  The  computer  interface  has  not  yet  been  constructed 
however.  If  the  computer  interface  is  developed,  the  system 
will  also  be  capable  of  measuring  multipoint  statistics  of 
any  high-speed  events  (up  to  100  MHz) . 

Experimental  measurements  have  been  performed  which 
verify  the  operation  of  the  photon  processor  with  100  MHz 
input  events  and  40  MHz  system  clock  operating  speed.  A 
sequential  photon  correlation  of  a low-level  sinusoidally 
modulated  optical  source  has  been  made  with  25  nsec  time 
resolution.  The  mean  velocity  of  a low-speed  jet  was  meas- 
ured with  the  dual  correlate  of  subtract  mode  under  conditions 
for  which  no  burst  counter  data  could  be  obtained. 

Photon  counting  experimental  work  requires  much  more  care 
in  selection  and  use  of  photomultiplier  tube  components  than 
does  burst-counter  work.  However,  single  photoelectron  pulse 
height  statistics  and  dynode  collection  efficiency  can  easily 
reduce  the  effective  signal  power  to  noise  power  ratio  by  a 
factor  of  4 for  two  tubes  with  the  same  quantum  efficiency. 

This  can  significantly  effect  a burst  counter  system  sensitivity. 

The  statistical  distribution  of  the  amplitudes  and  rates 
of  occurrence  of  classical  bursts  has  been  shown  to  be  central 
in  the  problem  of  specifying  or  predicting  the  data  rates  and 
errors  from  any  type  of  LV  signal  processor.  Differential  and 
cumulative  rate/amplitude  distributions  have  been  formulated 
and  analyzed  theoretically  and  have  been  measured  experimentally 
for  an  argon  backscatter  LV  system.  The  results  indicate  that, 
for  the  data  obtained,  the  smaller  aerosols  contribute  more 
to  the  phouon  correlation  accumulator  than  the  larger  ones. 

For  the  data  measured,  there  would  have  been  available  less 
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than  300  signals  per  second  adequate  in  magnitude  to  pro- 
duce burst  counter  data  from  scatterers  larger  than  0.7 
micron  in  diameter  while  there  would  have  been  over  100,000 
signals  per  second  producing  photon  resolved  signals  from 
0.2  - 0.3  micron  diameter  particles. 

We  have  demonstrated  both  theoretically  and  experi- 
mentally that  there  are  LV  signals  that  require  the  use  of 
photon  counting  techniques?  and  we  have  delivered  a new 
type  of  photon  counting  instrument  which  will  potentially 
operate  to  10  nsec  time  resolution  with  full  4 bit  multi- 
plication. 

Further  effort  with  computer  controlled  data  collec- 
tion is  required  to  explore  the  potential  of  the  new  system 
and  fully  define  its  limitations . 
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APPENDIX  A 

FILTERED  INHOMOGENEOUS  POISSON  PROCESSES 

In  this  appendix  we  have  provided  introductory  back- 
ground theory  for  filtered  inhomogeneous  Poisson  processes. 
For  more  introduction,  the  reader  is  referred  to  Papoulis 
[7]  or  to  the  more  advanced  text  by  Snyder  [63. 

A. 1 INHOMOGENEOUS  POISSON  IMPULSE  PROCESSES 

The  input  to  a random  linear  system  is  a inhomogeneous 
Poisson  impulse  process  z(t)  given  by 


z (t)  - l 6 (t  - t.)  (A.  1) 

— 00  A 

where  {t^}  is  the  set  of  random  occurrence  times,  X (t)  is 
the  instantaneous  statistical  mean  value  of  z(t),  (and  also 
the  mean  rate  of  occurrence  of  the  t^'s),  and  6(t)  is  the 
dirac  delta  function.  The  random  variables  are  independent 
of  each  other  statistically  and  obey  the  inhomogeneous  count- 
ing law,  i.e.,  the  probability  of  n * k occurrences  in  the 
interval  is 


P{n(tlft2) 


(A. 2) 


where 


y * / X (t) dt 
t. 


(A. 3) 
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The  quantity  y is  also  the  mean  and  variance  of  the  random 
variable  n(t^,t2). 

A. 2 THE  RESPONSE  OF  A RANDOM  LINEAR  SYSTEM- CAMPBELL ' S THEOREM 

The  output  of  the  random  linear  system  ■ (t)  * s the  super- 
position of  the  response  h(t-t^,  Y^)  to  each  input  impulses 

00 

a(t)  - l h (t  - x . , Y.)  (A. 4) 

— 00 

where  ( Y^}  is  a set  of  identically  distributed,  independent, 
vector  random  variables.  The  random  variable  Y affects  the 
shape  and  amplitude  of  the  response  function  h(t,Y).  In  the 
case  of  the  PMT  signal  it  may  take  the  form  of  a Bingle 
scalar  amplitude  variable.  In  the  case  of  the  classical 
optical  signal  from  turbulent  flow  both  a random  amplitude 
parameter  and  one  or  more  random  shape  parameters  due  to 
velocity  magnitude,  direction,  and  probe  volume  translational 
entrance  location  may  be  required.  The  theory  should  be 
applicable  so  long  as  *he  set  of  multidimensional  random 
variables  is  independent  of  the  set  of  occurrence  times 
{ t ^ • The  generalized  Campbell' s theorem  results  for  the 
instantaneous  statistical  mean,  variance,  and  auto-covariance 
of  s(t)  are  given  below,  they  apply  regardless  of  whether 
individual  pulseB  are  resolved  or  not. 

00 

<8  ( t)  > - / X(x)<h(t-T,  Y>dx  (A. 5) 

— oo 

00 

o2(t)  - <s2(t)>  - <s(t)>2  - / X(x)<h2(t- T,Y)>dx  (A. 6) 
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cov  s(t1)s(t2)  - <a (t^) a (t2) > - <8 (t^) ><s (t2) > (A. 7) 

00 

-J  X(T)<h(t1-  T,Y)h(t2-  T,Y)>dT 

“OD 

where  < > denotes  expectation  with  respect  to  *Y  inside  the 
integral  signs.  For  a causal  signal  such  as  that  from  the 
PMT,  where  hit)  is  zero  for  t < 0,  the  the  upper  limits  of 
integration  may  be  replaced  by  t or  the  minimum  of  t^  and 
t2  in  equation  (A. 7).  For  the  transient  case  where  the 
impulse  signal  z(t)  is  applied  at  t • 0,  the  lower  limits 
of  integration  may  be  replaced  by  0. 
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APPENDIX  B 


CAUSES  AND  EFFECTS  OF  THE  DISTRIBUTION  OF 
CLASSICAL  SIGNAL  BURST  AMPLITUDES 


The  second  phase  of  this  contract  was  an  attempt  to 
experimentally  and  theoretically  describe  the  statistics 
of  signal  burst  rate  and  amplitude  and  relate  these  to  the 
particle  size  distribution,  the  optical  parameters,  and  the 
electronic  parameters  in  a manner  which  would  ultimately 
allow  prediction  of  behavior  of  both  photon  counting  and 
burst  counter  data  processing  electronics. 

We  accomplished  our  objectives  for  some  specialized 
cases  of  small  particles  in  the  submic.ron  size  range  and 
somewhat  specialized  optical  geometry.  We  defined  rate 
and  amplitude  functions  which  should  be  studied  further,  and 
we  developed  the  theory  enough  to  conclude  that  with  data 
such  as  are  obtained  experimentally  (argon  backscatter)  the 
primary  contributors  to  photon  correlation  were  the  smaller 
amplitude  signals.  The  following  material  is  reproduced 
with  heading  and  reference  changes  from  our  Second  Interim 
Report  (an  unpublished  document)  without  further  refinement 
of  the  theory. 


B.l  RATE  AND  AMPLITUDE  DISTRIBUTIONS  OF  THE  CLASSICAL 
SIGNAL 


In  this  Appendix  we  are  concerned  with  the  probability 
density  for  the  signal  burst  amplitude  X j , the  rate  of  burst 
occurrence  R(t)  and  additional  rate/amplitude  distributions 
which  will  be  defined.  These  distributions  are  essential  in 
statistically  relating  an  accumulated  average  LV  measurement 
result  to  the  size  of  the  particles  which  contributed  to  the 
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result,  to  the  data  rates,  the  processor  efficiency,  bias 
errors,  and  other  quantities  of  interest.  We  are  not  aware 
of  any  previous  efforts  of  this  type. 

B . 1 . 1 The  Optical  Response  Function 

We  restrict  the  velocity  field  to  small  turbulence  with 
the  mean  flow  in  the  a direction  normal  to  the  fringes  of  a 

“a 

one-component  dual  scatter-system.  Off-axis  collecting  optics 
with  a rectangular  spatial  filter  are  placed  either  with  the 
optical  axis  in  the  plane  of  the  two  matched  transmitter  beams 
(case  1)  or  in  the  normal  plane  which  includes  the  bisector 
(z  axis)  of  the  transmitter  beams  (case  2).  In  either  case 
the  width  of  the  rectangular  aperture  in  the  image  plane  of 
the  collecting  lens  is  2cwq  where  c is  a constant  and  wQ  is 
the  1/e2  radius  of  the  transmitted  beams  in  the  probe  volume. 
The  projection  of  the  length  of  the  rectangular  aperture 
defines  the  length  of  the  probe  volume  in  the  z direction. 

For  case  ? *here  are  no  end  effects;  for  case  1,  scatterers 
passing  through  the  ends  of  the  probe  volume  will  produce  burst 
wave  forms  with  one  side  clipped  considerably.  The  length  Lz 
will  be  assumed  long  enough  to  neglect  these  end  effects,  but 
short  enough  to  maintain  high  fringe  contrast  and  uniformity 
of  intensity  in  the  z direction.  The  scatterers  are  assumed 
small  enough  for  validity  of  the  fringe  model  and  the  assump- 
tion of  a broad  scattering  pattern  so  that  perfect  signal 
fringe  contrast  is  obtained.  With  all  these  assumptions  we 
may  determine  the  optical  power  collected  at  the  PMT  photo- 
cathode from  a single  scatterer  at  position  (x ^ » y j , 2 j ) in 
terms  of  A as  follows 


,(x,y,z)  “ ye 


(1  + cos  2 tt  f^  x)  Rect  [ ] 


(B.l) 
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when 


2 | < Lz/2, and  is  zero  when  jz 


> L / 2 , where 

t* 


a = y , case  1 
a * x , case  2 

and  where 


2-0 
* ' ^ ? 


r • C&)  ff-)^ 


2 


fringe  period 

i(a)^  2P. 

— ?)  (— V) 

( 2irF j it w ^ 


(B.  2) 


D ■ collecting  lens  diameter 

T ■ combined  transmission  efficiency  of  transmitter 
and  receiver  optics  including  any  spectral  filter. 


aq  “ optical  wavelength 


F ■ focal  length  of  collecting  lens 
■ total  laser  power 

i(a)  ■ scattering  coefficient  (Van  deHulst)[21] 


The  expression  for  y is  broken  down  in  parentheses  to  sequen- 
tially show  the  effects  at  the  PMT  cathode,  the  collecting 
lens  area,  the  scattered  intensity,  and  the  incident  power 
density  at  the  center  of  the  probe  volume  (excluding  fringe 
variation  from  pedas^al.)  The  Rec.t  function  shows  the  effect 
of  the  width  of  the  spatial  filter  aperture.  The  Mie  scatter- 
ing coefficient  defined  by  Van  DeHulst  as  a function  of  a,  the 
particle  radius,  is  selected  as  i)(  or  ij.  depending  on  whether 
the  (linear)  polarization  of  the  incident  light  is  parallel  or 
perpendicular  the  plane  of  scattering. 


We  now  assume  that  the  jth  scatterer  passes  the  yz 

plane  at  time  t.,  and  at  location  (y-,z.)  with  velocity  Ua 
3 3 3 * 

and  scattering  coefficient  i (a^ ) . This  results  in  the  time 

response 
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X . (t)  ■ X.  Rect[a/2cw  ] Rectlz./L  ] 

J J ® J Z 

-2  [ (t-T  • )U]2/w2 

e -1  °[1  ♦ cos  (2nf  (t  - t . )U] 

X J 


(B.3) 


a « y. 


case  1 


where 


(t  - t j )U,  case  2 


-2yj/w2 
■ v(a.)e  -1 

B.1.2  Probability  Density  of  X. 


(B-4) 


The  result  of  equation  B.3  is  the  product  of  a normalized 
function  displaced  in  time  to  t - Tj  with  a random  amplitude 
which  is  the  product  of  two  independent  random  variables 
Yj  * t C a j ) and  Yj  where 


Yj  “ fl(aj} 

Y j - fjCYj)  " exP  l -Zy^/w^lRect  [^-]Rect  [j-i] 


(B.  5) 


In  a more  general  model,  the  signal  amplitude  would  also 
depend  on  direction  and  other  position  variables  (x.,y.)  and 

J J 

the  shape  of  the  normalized  function  would  also  depend  on 
thesj  quantities.  The  purpose  of  our  simplified  model  is  to 
examine  the  effect  on  px  (x)  of  the  uepender.e  of  x^  on  both 
particle  radius  and  probe  volume  entrance  location.  Papoulis 
[7]  gives  theory  required  to  obtain  density  p,(x)  for  the 
product  y j y j ■ Xj  in  terms  of  the  density  p^M  and  Py(Y) 


PX(X)  " L ITT  PT^)PY(a)da 

" -L  TTf  PY(a)PY(a)da 


(B. 6a) 


(B. 6b) 
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where  in  B. 6 a is  a dummy  integration  variable.  The  densi- 
ties p (y)  and  Py  are  straight-forward  in  principle  to  obtain 
from  the  particle  radius  a^  density  and  the  density  for  the 
location  y^ . Papoulis  gives  the  density  for  a function  of  a 
random  variable  from  which  we  obtain: 


P„(y)  - 


PQ  (a) 


at  a ■ f j (y) 


(B.7a) 


Py(Y) 


pv(y) 


at  y ■ f2  1 (Y) 


(B  - 7b) 


The  notation  in  these  formulas  means,  for  example  in  B.7a, 
the  argument  of  the  density  pa (a)  is  replaced  by  the  inverse 
function  of  y.  This  density  is  divided  by  the  first  derivative 
of  the  function  f^  (where  y * f^(a))  and  evaluated  at  a ■ in- 
verse function  of  y* 

When  we  consider  the  evaluation  of  equations  B.6  and  B.7 
even  for  our  simplified  model  given  by  equations  B.l  through 
B.5,  we  see  that  the  relationship  of  signal  amplitude  density 
to  particle  diameter  density  is  very  complex  and  is  not  related 
by  a one-to-one  function  because  of  the  random  entrance  loca- 
tion multiplier  Y.  Complete  evaluation  would  require  the  Mie 
scattering  function,  its  first  derivative,  its  inverse  and 
the  particle  size  density. 

B . 1 . 3 Simplified  Evaluation  Examples 

The  nature  of  the  results  in  equations  B.6  and  B.7  can 
be  demonstrated  by  the  use  of  simplified  examples.  First  we 
consider  monosized  scatterers  and  then  we  consider  a rth 
power  law  approximation  of  the  Mie  scattering  function. 
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B. 1.3.1  Monosized  Scatterers 

With  a ~ a for  all  scatterers,  we  obtain 
o 


Y ■ Y, 


Py(Y)  a 5(y  ' Y0) 


1 


where  6 is  the  dirac  delta  function.  This  makes  evaluation 
of  equation  B.6  simple  once  we  determine  Py(Y).  First  con- 
sider case  1 where  the  receiver  aperture  restricts  |y^|  to 
values  less  than  cwQ.  The  function  Y(y)  is  even  and  we 
obtain  correct  results  by  considering  the  half  interval 
0 < y.  < cwo.  For  uniform  particle  distribution  in  volume 
with  no  velocity  gradients,  the  density  for  y.  is  uniform 
and  given  by 


Pyi(y) 


* . 0 < y < cw„ 

cwQ  7 o 

■ 0,  otherwise. 


From  B.  7a  and  B.  5 


Py(Y) 


2/1  cY/^InY 


- 2c* 

, e < Y < 1 


otherwise 


This  density  is  plotted  in  Figure  B1  with  c - 1 (aperture 
matched  to  1/e  probe  width.)  The  dashed  line  indicates 
the  shape  (values  off  by  1//7)  which  would  result  with  the 
aperture  matched  to  the  1/e  radius  wq//7. 
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When  equations  B.g  and  B.10  are  substituted  into  B.6a 


we  obtain  the  density  p.(x)  as 


p*(x)  - r py  (x/V 


(B. 11) 


2 /?  c X /- ln(x/YQ5 


-2c 

, ye  < X < y, 


In  other  words,  with  monosized  particles,  the  density  for 
the  burst  amplitudes  X^  assumes  the  shape  of  the  density  for 
the  random  probe  volume  multiplier  Y.  When  the  case  2 
optical  system  is  considered,  the  result  is  the  same  except 
c is  chosen  much  larger  than  1.  When  this  is  done  the 
results  become  singular  in  the  limit  of  large  c.  Qualitatively, 
the  result  shows  large  relative  probability  of  x^  near  zero. 

We  will  show  later  how  the  singularity  vanishes  in  application 
problems . 


B . 1.3.2  Power  Law  Scattering 

When  the  distribution  of  particle  diameter  lies  entirely 
in  the  Rayleigh  size  regime,  then  the  Mie  scattering  (equation 
B.2)  behaves  as  the  6th  power  of  the  radius  a.  In  some 
cases  we  may  approximate  the  function  y(a)  as  an  rth 
power  function  over  a significant  range  of  values  of  a so 
that 


Y (a) 


o 


(B  . 1 2 ) 


where  yq  is  the  value  of  y at  a selected  value  of  a ■ ap. 


Let  us  assume  that  in  the  same  range  of  a,  the  density  p (a) 

3 


for  particle  radius  obeys  a qth  power  law: 

q 


?aM  ' Pao(4’ 


(B.131 


1 i 


i i 
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l 

where  paQ  is  a constant.  In  this  case  we  may  evaluate  equa- 
tion B.7a  to  determine  the  behavior  of  py(y)  in  this  range. 

The  result  is  i 


(l+q-r)/r  ; 

p (Y)  - p f-I-)  CB- 14) 

‘yo'-y0  | 

4 

1 

where  py£)  is  a constant.  As  an  example,  if  Rayleigh  • 

scatterers  (t  ■ 6)  were  distributed  with  a density  behaving 

- 1 • S 1 

as  q - -4,  then  the  density  for  y behaves  as  (y/y0)  . i 

^ j 

By  comparison,  with  r ■ 4 and  q ■ -7,  the  density  for  y would  ^ 

go  as  (y/y0) "2‘ S.  j 

In  order  to  determine  the  shape  density  pA(x),  given  \ 

Py(Y)  over  a wide  of  y,  we  could  again  utilize  equation  B.6. 

This  has  not  been  done  at  this  time.  However,  the  results  of 
the  preceding  example  show  that  a particle  with  radius  a con- 
tributes a value  of  X between  e"2c  Y(a)  and  y(®).  The  results 
of  the  present  example  are  referred  to  in  later  sections. 

B.1.4  Burst  Rate 

B.l.4.1  Introduction.  The  nature  of  an  LV  signal  is 
affected  by  all  of  the  burst  signals  present  regardless  of 
whether  they  are  of  sufficient  amplitude  to  appear  as  classical 
signals  and/or  cross  a threshold  level.  Theoretically  we  can 
work  with  X(t)  even  in  the  photon-resolved  cases,  and  must 
therefore  have  precise  definition  of  the  rate  R(t)  of  occur- 
rence of  the  bursts  in  addition  to  the  probability  density  ! 

p (x)  for  their  amplitudes.  If  we  were  given  monosized 
scatterers  with  number  density  p (scatterers/m'5)  and  a sharply 
defined  probe  volume  cross-sectional  area  Ay  normal  to  the 
vector  velocity  V then  obviously 

R(t)  - | VCt)  | p Ay  (B.15)  i 
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If  velocity  gradient  effects  were  to  be  included,  a surface 
integral  would  replace  equation  B.1S,  and  the  formulation  of 
A(t)  as  an  inhomogeneous  filtered  Poisson  process  would  have 
to  be  extended  as  a superposition  of  such  processes. 

In  many  practical  problems,  some  of  the  boundaries  of 
the  probe-volume  cross-sectional  area  are  not  sharply  defined 
and  an  arbitrary  definition  (such  as  the  1/e  intensity  con- 
tour) is  meaningless  when  a wide  distribution  of  values  of  y 
(particle  size)  is  present.  In  such  cases  equation  B.15  must 
be  used  with  definition  of  p and  Ay.  The  cumulative  number 
density  of  particles  with  radius  a greater  than  a'  is 


N_ (a  > a' ) ■ / “ n (a)da 


(B  .16) 


and  p • N (a>o)  where  n (a)  is  a differential  size/number 

B cl 

density  distribution.  We  observe  that  p fa)  andn_(a)  are 

B B 

related  by  the  normalizing  constant  p 


PaCft) 


Ulil 

p 


(B  .17) 


Thus,  the  relationship  between  the  cumulative  probability 
function  and  the  cumulative  number  density  function  is 


p(a<a' ) ■ 1 - jNa(a>a') 


(B.18) 


B.  1.4.2  Burst  Rate/Amplitude  Distribution.  When p and 
Ay  are  well-defined  then  equation  B.15  is  applicable.  We 
may  also  define  a differential  distribution  r^(t,X)  given  by 


rA  (t, A)  ■ R(t)  pA  (x) 


(B.19) 
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The  cumulative  rate  is  the  integral  of  this  function: 


Rx (t , X > X * ) 


rx(t,x)dx 


(B. 20) 


\A 


in  this  definition,  R(t)  ■ R^(t,x>o).  The  definitions  of 
the  rate  functions  seem  redundant.  However  with  extended 
definition  as  a limit  of  a conditional  function  (Sections 
B.  1.4.3  and  equation  B. 32) , there  are  cases  where  rx(t,X) 
has  meaning  where  p x ( X ) does  not. 

B.l.4.3  Extended  Definitions  of  Number  Density  and  Cross 
Sectional  Area 

In  practice  we  must  deal  with  experimental  situations 
where  it  is  not  practical  to  determine  the  shape  of  Px(x) 
precisely  in  the  range  of  small  values  of  X^.  In  some  such 
cases,  the  signal  and  noise  contributions  from  the  small  sig- 
nals  may  be  lumped  together  approximately  as  an  additional 
background- light  term  without  significant  error,  but  this  is 
not  always  true.  An  approach  is  needed  which  allows  the 
analysis  to  proceed  to  the  point  of  ascertaining  the  effects 
contributed  by  the  small  X tails  of  the  density  p^(X).  We 
find  this  in  the  use  of  conditional  densities  and  rates.  We 
assume  that  all  particle  radii  "a"  are  greater  than  some 
minimum  radius  a^  and  that  the  probe-volume  cross-sectional 
area  is  limited  by  sharply  defined  boundaries  in  terms  of  a 
parameter  ft,  ("c"  in  the  preceding  optical  examples) . Then 
we  determine  p C A | a ^ , S’) , use  it  in  the  analysis  of  the 
processor  behavior  to  obtain  a result,  and  finally  take  limits 
as  a^  goes  to  zero  and  the  cross  sectional  area  goes  to 
infinity.  In  the  last  step,  precise  models  for  pfl(a)  and  the 
optical  response  function  only  become  necessary  if  the  final 
result  shows  significant  changes  as  the  limits  are  evaluated. 
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In  our  case  2 optical  example  we  implicitly  demonstrated 
the  approach  of  finding  pA(X)  conditioned  on  c.  In  order  to 
find  the  conditional  density,  given  that  all  scatterers  are 
larger  in  diameter  than  a, , we  replace  p (a)  in  the  formulas 
with  pa(a|a>a1)  given  by 

n (a) 

Pa(./a  » »j)  - nt;,>ti)-  . *>»! 

■ 0 a <_  a^ 

instead  of  using  (B.17).  Similarly,  p ■ N_(a>o)  is  replaced 
by  Nfl(a>a^)  in  all  formulas.  For  example,  equations  B.19 
and  B.15  become 

rA(t,X,a1,X)  - R(t|a1,X)pX:j  (X  |a:,  £)  (B.22) 

R(t|alfK)  - |V(t)|  Na(a>ai)Av(X)  (B.23) 

when  the  limit  of  r A (t , X ,a^  ,X)  is  determined  as  a^  0 and 
Ay  •*  « as  in  equation  4.9,  then  the  function  rA(t,x) 
which  results  is  well  behaved  except  possibly  at  X ■ 0,  and 
RA(t,X  > x')  is  defined  for  all  values  of  X*  > 0.  Physical 
reasoning  shows  that  the  product  XR(X>X')  approaches  a 
finite  or  zero  value  as  X ■+  0,  since  this  product  represents 
the  signal  power  contained  in  the  low -ampl itude  bursts. 
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B.2  APPLICATION  TO  PHOTON  CORRELATION 

. We  are  now  able  to  quantitatively  describe  the  expected 

* value  of  the  output  of  the  accumulators  in  a photon  correla- 

- tion  experiment  in  terms  of  the  system  parameters.  We  also 

have  a way  of  looking  at  how  the  accumulator  output  is 
f weighted  with  respect  to  the  distributions  of  the  burst  ampli- 

f tudes. 

t 

B . 2 . 1 Expected  Accumulator  Values 

1 

We  may  rewrite  equation  2.6*  in  conditional  notation 

with 


<nk|A(t)>  * At  A (kAx)  (B.24) 

Following  through  with  A(t)  • nP(t)/hv,  the  expected  value 
of  the  accumulator  sum  is 

N 

<Aj>  ■ ‘Jj’Wj  ' N<nknk*j>  (B-25) 

* N(iT)2<X(t)X[t*jiT)> 

when  j + 0.  Equation  (2.15)*  omitted  background  light  and 
variation  of  the  shape  of  the  envelope  of  f(t)  with  w^. 

These  effects  are  now  included  in  equations  3.1  -3.8,  when 
the  rate  of  occurrence  R(t)  of  the  bursts  is  a constant  R. 

For  the  case  1 optical  model  of  equation  B.3  and  B.4 

a — 

In  this  subsection,  an  asterisk  on  an  equation  number  ref,  . s 
the  reader  to  Interim  Report  #1,  August,  1974.  That 
report  is  unpublished,  but  the  material  is  not  necessary  here. 


113 


AEOC-TR-76-81 


with  steady  x directed  flow  we  obtain: 


(B. 26) 


R<xj>/7  W0  .,V/«*  , 


ll  + J cos ( 2 it f xUt ) ] 


When  we  recall  that  for  V » Uax  then  the  burst  rate  parameter 
becomes  from  equation  B.1S 


R - UpA, 


(B. 27) 


where  p was  the  scatterer  number  density  and  Av  is  the  crossec 
tional  area  used  to  find  px(x).  Thus  we  see  that  while  the 
velocity  U remains  in  the  shape  of  <X (t) X (t  + x) > , it  cancels 
in  the  amplitude  multipliers.  If  we  note  that  in  B.25 


(B. 28) 


where  T is  the  total  data  collection  time,  the  result  of 
combining  equation  B.25  - B.28  is 

<Aj  > • TAt(Xb  +<Xj>pAv  WQ/j  ) 

<X?>TAtpA  w /?  -t7'U2/w„  j 

+ — i— j — —Z—  e °ll  + J cos(2itfxU  t)] 

This  result  is  related  to  fundamental  system  parameters 
through  p (X)  by  the  equations 


*This  assumes  that  the  pedastal  and  AC  signal  are  non-over- 
lapping spectrally  so  that  cross-correlation  terms  between 
pedastal  and  AC  portions  vanish.  It  also  neglects  the  behav- 
ior in  the  vicinity  of  x ■ o (x^  <<  At). 
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<xj  * " ) „ >px (X)dA  (B.30) 

and 

<*?>  - / *x2p  (x)dx  (B . 31 ) 

J . CB 

B.  2.2  Source  Distribution 

In  equation  B.  29  we  used  the  simple  definition  of  burst 
rate  given  by  equation  3. IS  and  avoided  the  complexities  of 
Section  B.l.3.3.  Generally  speaking,  the  function  p x C x ) has 
no  useful  definition  independent  of  Ay;  likewise,  p is  an 
ambiguous  term  because  there  may  be  large  numbers  of  scatterers 
too  small  to  contribute  significantly  to  the  accumulator  out- 
put. The  significant  quantity  in  B.29  is  the  limit  function 
C(X)  defined  by 

r.  ( x, a,  ,K) 

«(*)  • lim  — — (B.  32) 

A,,  -*■  “>  u 


where  equation  B.22  and  B.23  are  used  with  |V(t)|  ■ U.  When 
equation  B. 32  is  substituted  in  B.29  with  the  identification  of 
W ■ w /tt/2  as  the  effective  probe  width  in  the  direction  of 
flow  we  obtain 


W o® 

<A . > - TA t [ (X.  + W / " xc(x)dx)2  + — / X2e(X)dX  • shape  ) (B.33) 
J 0 o /7  o function 

Equation  B.33  applies  to  the  result  of  using  a real-time 
digital  correlator  with  full  multiplication.  In  the  dual- 
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correlate  and  subtract  system  we  are  to  construct,  only  a 
portion  of  the  right  hand  integral  form  in  B.33  remains  in 
the  accumulator  result;  but  the  constant  term  contributes 
to  the  statistical  variance  of  the  result.  The  quantity 
which  determine  the  weighting  of  the  result  according  to 
the  amplitudes  X^  and  their  relative  number  is  X C(X). 

This  weighting  function  i9  well  behaved  even  if  e(X)  has  a 
first  or  second  order  pole  at  X ■ 0,  and  if  t(X)  becomes 
monotonically  decreasing  faster  than  x beyond  some  value  of 
x. 

If  the  product  X £(X)  is  still  significant  at  values  of 
x where  single  photoelectron  pulse  overlap  or  discriminator 
dead  time  effects  begin  to  dominate,  then  these  effects  will 
truncate  the  weighting  factor  for  the  larger  volumes  of  X. 

If  5(i)  is  not  decreasing  faster  than  X at  this  point  the 
maximum  of  the  truncated  weighting  function  will  occur  at  the 
point  where  nonlinear  saturation  effects  set  in.  We  have  no 
idea  at  this  time  what  the  effect  of  this  region  will  be  on 
the  output. 

B . 2 . 3 Example;  Monosized  Scatterers 

The  use  of  the  notation  ?(X)  simplifies  notation  while 
obscuring  the  communication.  As  an  enlightening  example  we 
return  to  the  case  2 optical  system  with  monosized  scatterers. 
In  that  case  the  use  of  p for  number  density  was  not  ambiguous, 
and  Av  can  be  conditionally  defined  as 

Av  ■ 2cwqLz  . (B.34) 

when  this  is  done,  equation  ti . 1 1 gives  the  conditional  density 
p^  (X  | c)  and  we  have 
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lim  p(2cwQL2)p 


(i  |c) 


(B.35) 


W0LZ  p * 


0 < x < Y, 


where  yQ  is  the  peak  value  of  x^  for  particles  through  the 

center  of  the  probe  volume.  The  definition  of  $(x)  avoids 

the  anomalies  which  arise  in  the  definition  of  p.  (X)  due  to 

A 2 

probe  volume  boundaries  which  are  not  well  defined.  The  x 
multiplier  (or  the  X multiplier  in  determining  the  constant 
part  of  <Aj>)  removes  the  singularity  at  x ■ 0.  The  function 
X//-ln(X/Y^)  is  plotted  in  Figure  B2  foT  yQ  ■ 1. 

B.2.4  Interpretation  for  Particle  Size  Distribution 

If  we  recall  equation  .6a  and  insert  this  in  the  defini- 

2 2 

tion  of  the  weighting  factor  x C(X)  we  see  that  the  X term 
may  be  carried  inside  the  integral  with  the  result: 

X25(X)«  lim  / PY(«)l  — r^T"  pY(«  |Av)lda  (B.  36) 

A -►«  - » T 'I 


We  recognize  this  as  a superposition  integral  with  the  term 
in  brackets  being  the  weighting  function  for  a single  particle 
size  which  produces  a value  of  y ■ a,  The  integral  adds  these 
functions  with  weighting  by  the  density  for  y. 

We  have  finally  reached  the  point  where  we  can  determine 
the  composition  of  the  photon  correlator  accumulator  result 
in  terms  of  particle  size.  To  do  this  we  reinsert  X25(X)  from 
equation  B.36  into  equation  B.33  and  interchange  the  order 
of  integration.  Integrating  with  respect  to  X first  leaves 
the  second  term  in  B.33  as 
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W a,  00  pA  X . 

— ! P (a)  [ lirn  / p (A|A  )dx]da  CB.  37) 

n 7 Av-»“  o 11  v 

2 

A similar  expression  with  X instead  of  x holds  for  the 
first  term  in  B.33.  The  result  of  the  inner  integral  in  equa- 
tion B.37  is  a function  of  the  dummy  integration  variable  a 
which  represents  y . The  variable  y is  directly  related  to 
particle  size  through  the  Mie  scattering  law. 

The  theory  which  relates  the  particle  size  distribution 
to  the  photon  correlator  accumulator  outputs  is  more  complex 
than  we  originally  believed.  We  have  obtained  results  theo- 
retically, but  at  the  present  we  have  not  had  time  to  extend 
the  generality  beyond  a specific  somewhat  simplified  optical 
model  or  to  evaluate  the  integrals  parametrically  to  obtain 
more  graphic  results. 

B.3  APPLICATION  TO  BURST  COUNTER  SYSTEM 

The  objectives  of  the  first  extension  of  our  present 
contract  included  development  of  a statistical  signal  model 
that  would  benefit  both  the  development  of  photon  correlation 
and  the  further  refinement  of  burst-counter  systems.  We  can 
not  go  much  beyond  opening  the  door  with  inspect  to  burst  - 
counter  applications  because  of  scope  limitations;  but  we 
would  like  to  briefly  indicate  the  nature  of  several  tentative 
conclusions  we  have  drawn  without  having  had  time  to  pursue 
very  detailed  analytical  or  computational  studies. 

B . 3 . l Transient  Filter  Response-SNR 

Wo  have  experimentally  observed  filtered  low- level  IV 
signals  which  not  only  demonstrated  the  non-stationary  and 
signal-dependent  nature  of  the  noise  but  also  exhibited  a 
time-delay  shift  of  the  high-frequency  portions  of  the  noise 
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with  respect  to  the  signal  waveform.  In  this  observation, 
the  positive-going  zeio-crossing  portions  of  the  signal 
apperred  to  be  much  less  noisy  than  the  negative-going  zero- 
crossing portions.  In  this  instance  it  was  significant  that 
the  signal  polarity  was  chosen  to  count  fringe  periods  on 
the  positive- going  zero-crossings  and  not  the  negative-going 
zero-crossings . 

When  Gaussian  noise  is  filtered  by  a linear  filter,  only 
the  magnitude  of  the  frequency  transfer  function  for  the  fil- 
ter is  significant  for  the  output  noise  of  the  filter.  This 
is  true  because  the  statistics  of  the  output  noise  are  com- 
pletely determined  by  the  output  power  spectrum  which  is  the 
product  of  the  input  power  spectrum  and  the  magnitude  square 
of  the  Fourier  transform  of  the  filter  impulse  response. 

The  experimental  observation  mentioned  in  the  above  paragraph 
and  the  theory  of  Section  2 show  that  the  time  distribution 
of  the  noise  with  respect  to  the  signal  is  affected  by  the 
transient  response  of  the  filter  (or  its  entire  Fourier 
description,  including  the  phase  characteristic.)  This 
matter  should  be  pursued  further  in  the  future. 

B.3.2  Detection  Improvement 

In  a burst-counter  the  functions  of  detecting  the  pres- 
ence of  a burst,  estimating  the  mean  period,  and  dec iding 
whether  the  data  word  is  acceptable  are  three  separate  func- 
tions. With  state-of-the  art  error  rejection  (decision) 
circuits,  valid  data  can  be  obtained  by  setting  the  detection 
threshold  low  and  rejecting  invalidated  data.  However,  informa 
tion  is  lost  where  the  threshold  detection  circuit  triggers 
on  a noise  fluctuation  and  the  estimation  dead  time  causes 
the  system  to  miss  a burst  which  would  have  produced  a vali- 
dated measurement. 
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There  are  two  simple  ways  to  significantly  improve  the 
threshold  detection  efficiency  of  a burst-counter.  The  first 
is  to  use  state-of-the  art  PMT's  to  minimize  probability  of 
occurrences  of  anomalously  large  single-photo-electron  pulses. 
The  second  is  to  provide  a separate  signal  path  and  filter 
for  the  threshold  detection  and  the  period  estimation  circuits. 
To  be  specific,  if  the  optical  system  restricts  the  field  of 
view  of  the  PMT  so  that  all  or  almost  all  pedastal  signals 
are  accompanied  by  a sinusoidal  signal-burst  component,  then 
a good  detection  filter  is  a low-pass  pedastal  detection 
filter  with  much  less  bandwidth  than  is  required  for  the 
period  estimation  circuits. 

B.3.3  Particle  Size  Effects 

For  a given  LV  system  and  set  of  measurement  parameters 
we  have  the  following  composition  for  the  output  validated 
data  rate: 


Rvd  ■ /*  Vx>V(X,lU  * Pv.J  <B-38’ 

where  r^(X)  is  the  burst  rate  amplitude  distribution  described 
in  Section  B.1,  Pgd(X)  is  the  probability  that  a burst  of  ampli- 
tude X.  - x will  be  detected  and  produce  a validated  estimate, 
R , is  the  rate  at  which  the  processor  triggers  on  background 

£U 

fluctuations  (estimation  produces  erroneous  and  irrelevant 
data),  and  P d is  the  probability  that  validation  occurs  with 
a false  trigger.  The  quantities  Pgd,  Red,  and  Pved  are 
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functions  of  the  background  light  level,  the  detection 
threshold,  the  processor  dead  time  and  the  total  burst  rate 
(things  get  much  more  complicated  when  the  probability  is 
appreciable  that  the  detection  threshold  will  be  crossed 
more  than  once  during  the  processor  dead  time;  also  if  there 
are  subthreshold  bursts  present  during  a detected  burst, 
the  noise  and  interference  are  higher  and  Pgj(X)  is  lower.) 
Generally,  the  problem  we  wish  to  solve  is  the  maximization 
of  the  good  data  rate  with  constraints  on  R .,  P , and  with 
other  constraints  (such  as  no  seeding,  a certain  system  cost, 
a minimum  acceptable  data  rate,  etc.)  The  statistical  signal 
model  we  have  presented  above  with  generalization  of  the 
optical  model  and  proper  modeling  of  the  flow/particle  dynamics 
and  the  processor  components  form  the  basis  for  constructing 
a system  optimization  code  which  could  be  used  in  measurement 
planning  or  on-line  system  optimization. 

The  complexity  of  equation  B.38  is  considerable  when 
all  the  dependence  of  the  many  deterministic  and  random  sys- 
tem parameters  are  shown  explicitly.  Some  of  the  contract 
time  was  devoted  to  attempts  at  being  more  specific  and  obtain- 
ing a very  crude  model  for  P ( X ) so  that  there  would  be  some 
quantitative  comparison  between  the  source  distribution  for 
photon  correlation  output  and  burst  counter  output.  The  task 
was  quite  beyond  the  limitations  of  time  under  which  we  were 
working.  The  significant  thing  we  have  done  in  this  area  Is 
establish  a common  approach  for  the  comparison  and  optimiza- 
tion of  LV  processors.  The  basis  is  in  the  function  r^(X) 

(or  the  related  £(X))  and  the  composition  of  this  function 
in  terms  of  the  particle  size  distribution  and  the  random 
probe  volume  multiplier  effects. 
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B.  4 MEASUREMENT  AND  SCALING  OF  RATE /AMPLITUDE  FUNCTIONS 

We  have  shown  that  measurement  time  (or  data  rates) 
for  photon  correlation  (or  burst  counters)  depend  on  an 
integral  superposition  of  £(X)  (or  r^(X))  which  has  been 
related  to  particle  size  number  density  and  optical  param- 
eters. In  this  section  we  consider  the  two  separate  prob- 
lems of  experimental  measurement  of  this  function  and  its 
scaling  behavior  with  changes  in  system  parameters.  These 
two  problems  are  related  because  the  quantity  x(t)  and  the 
separate  burst  amplitudes  Xj  are  not  always  directly  meas- 
urable, and  scaling  may  be  essential. 


B . 4 . 1 Direct  Measurement 

The  measurement  of  X(t)  exactly  is  not  possible.  The 
trade-off  is,  as  usual,  between  resolution  (bias  error)  and 
noise  (statistical  variation).  Equation  2.7  may  be  rewritten 


<i(t)>  = e <g£>  X(t) 


(B . 39) 


when  h(t)  is  a low-pass  pulse  and,  t ^ is  small  enough  for 
h(t)  to  behave  as  an  impulse  function  with  respect  to  the 
time  variation  of  X(t).  Under  the  same  condition  equation 
2.8  becomes  approximately 


°i2  ^ <gi.2> 

1 Th  1 


x (t) 


(B . 40) 


and  the  instantaneous  fractional  rms  deviation  is 


<4(1) 


<gi> 


(B.41) 
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Thus,  if  X ( t ) t >>  1,  i(t)  is  proportional  to  X(t)  with 
small  variation.  If  this  condition  were  valid  for  the  LV 
signals  of  interest  we  wouldn't  be  concerned  with  any  of 
the  present  effort. 

Fortunately,  for  "sparse"  bursts,*  the  requirements 
for  measurement  of  the  distribution  of  burst  amplitudes  X^ 
are  not  as  stringent  as  those  for  continuous  estimation  of 
X(t) . We  will  only  consider  a simple  sub-optimal  approach 
since  the  optimal  measurement  of  5X(X)  could  be  the  subject 
of  a study  in  itself.  If  a low-pass  filter  is  used  which 
passes  the  burst  pedastal  but  not  the  sinusoidal  signal, 
then  the  value  of  in  B.40  is  increased,  the  convolu- 
tion in  2,8  removes  the  sinusoidal  portion  of  x(t),  and  we 
would  replace  x(t)  in  R.41  with  Xp(t),  the  pedestals  only. 
From  equation  B . 3 we  note  that  this  leaves  X^  as  the  peak 
value  of  a Gaussian  pulse  at  the  point  where  the  fractional 
rms  deviation  R.41  is  minimum.  Due  to  the  increase  in  th 
obtained  by  reducing  the  filter  bandwidth,  smaller  values 
of  X.  can  be  measured  than  if  a peak  detector  were  applied 
to  the  direct  output  of  the  photomultiplier  tube. 

In  order  to  extend  the  measurement  to  smaller  values  of 
Xj  there  are  several  straight-forward  improvements  of  the 

*The  phrase  "an  average  of  less  than  one  scatterer  in  the 
probe  volume"  is  very  ambiguous  here.  What  is  usually 
implied  in  the  LV  literature  is  that  for  the  detection 
threshold  and  filter  which  are  being  used,  the  average  burst 
separation  is  much  greater  than  the  burst  duration.  A typi- 
cal probe  volume  for  high  velocity  air  flow  measurements 
can  easily  contain  an  average  of  more  than  one  scatterer  if 
all  the  smaller  scatterors  are  counted.  With  typical  filter 
bandwidths  for  transonic  measurements,  the  bursts  produced 
by  the  smaller  scatterers  are  photon  resolved;  thus  the 
appearance  of  the  signal  is  that  of  "sparse"  bursts  plus 
background  when  in  fact  some  of  the  "background"  is  low- 
level  overlapping  or  non  over-lapping  signals. 
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above  procedure.  The  first  is  to  reduce  the  aiT  velocity 
to  a very  small  value  and  thereby  reduce  the  required  band- 
width and  extend  In  this  manner  values  of  Xj  may  be 

measured  by  a peak  detector  even  though  the  same  scatterer 
size  will  produce  photon-resolved  signals  at  higher  veloci- 
ties (with  wider  bandwidth  filters).  We  have  used  this 
method  in  the  experimental  work  reported  in  the  section  5.0. 

For  on-line  measurements  at  high  velocity  or  for  situa- 
tions where  the  particle  number  density  is  too  large  for 
the  "sparse”  condition  to  be  valid  another  approach  would 
be  to  expand  the  laser  beams  incident  on  the  transmitting 
lens  and  reduce  the  collecting  system  spatial  filter  aper- 
ture. This  would  increase  the  incident  optical  power 
density  by  the  square  of  the  expansion  while  increasing  the 
filter  bandwidth  required  only  linearly.  The  scaling  laws 
of  the  next  section  could  then  be  applied. 

Additional  improvement  of  the  above  approaches  would 
result  from  reducing  the  bandwidth  further  by  replacing  the 
filter  with  an  integrator  and  logic  which  measured  the  step 
level  change  (total  burst  energy)  and  then  reset  the  inte- 
grator to  zero  before  the  next  burst  occurred.  An  inte- 
grating particle  size  instrument  has  recently  been  developed 
by  M.  Farmer, J.  Hornkold,  et.  al.  at  SAI  Tullahoma.  A modi- 
fication of  that  instrument  or  similar  integrating  circuitry 
could  be  useful  in  later  efforts  along  these  lines. 

All  of  the  above  methods  become  unreliable  at  levels  so 
low  that  background  fluctuations  are  detected  as  signal 
bursts,  and  further  study  of  the  statistical  estimation  of 
5^(X)  by  photon  counting  methods  may  later  be  appropriate. 

B .4 . 2 Scaling 

The  scaling  laws  are  implicit  in  the  theory,  but  are 
stated  here  for  convenience.  With  the  exception  of  velocity 
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effects  (since  rx(X)  ■ |V|cx(x))  the  functions  rx(X),  and 
CA(X))  scale  the  same  way  so  we  refer  only  to  rx(x)  and  its 
integral  RX(A  > X^). 

B . 4 . 2 . 1 Velocity 

CA(X)  is  not  affected  but  our  ability  to  measure  the 
function  is  affected.  rx(x)  and  Rx ( x > x^)  are  proportional 
to  velocity  unless  velocity  change  also  changes  the  size/ 
number  density  distribution  of  scatters  in  the  air  (such  as 
the  occurrence  of  condensation  during  a wind  tunnel  test). 

B.4.2.2  Effective  Collected  Power 

An  increase  in  collected  effective  power  (laser  power, 
PMT  quantum-collection  efficiency,  optical  transmission 
efficiency,  collecting*  f/no)  by  a factor  k will  produce 
r A (x)  given  by 


r^X)  - £ rx(x/k)  (B. 42) 

The  1/k  multiplier  does  not  enter  in  the  integral  quanti- 
ties , i . e . : 

rx(X  > Xi)  • Rx(\  >-£=•)  (B.  43) 

This  is  important  since  the  cumulative  function  may  be  meas- 
ured and  differentiated  to  obtain  the  differential  function. 

B . 4 . 2 . 3 Incident  Beam  Diameter 

An  increase  in  the  waist  diameter  of  the  parallei.be'ams 
incident  on  the  transmitting  lens  by  a factor^of'^k  with 
corresponding  reduction  by  a factor  k "file  probe  diameter 


^Simplistic  model  foy -eirtafl  scatterers  not  considering 
scatter  lobe  patterns. 
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w , and  all  dimensions  of  the  collecting  system  spatial 
filter  aperture,  changes  the  shape  (number  of  fringes  of 
the  burst)  and  effects  the  distribution  functions  as 

rx(X)  " “T  rX 

R'x(x  > xl}  " “7  R(x  >^) 

In  this  instance  the  distribution  expands  to  the  higher 

values  of  X as  k because  of  the  square  increase  in  incident 

2 

intensity.  There  is  an  additional  1/k  factor  due  to  the 
square  decrease  in  probe-volume  cross-section  area.  We 
should  also  point  out  that  for  a photon  correlation  system 
this  change  also  reduce*  the  value  of  Wx  in  equation  B.33 
by  1/k. 

Comment:  If  the  collection  optics  spatial  filter  aper- 

ture is  not  scaled  along  with  the  probe  volume  beam  diameter, 
the  shape  of  rx(x)  is  changed  so  that  no  simple  scaling  is 
possible . 

B . 4 . 2 . 4 Transmitter  Focal  Length 

If  the  transmitter  focal  length  is  reduced  by  a factor 
k,  and  the  collecting  optics  spatial  filter  is  reduced  in 
uach  dimension  by  k,  and  the  scattering  angle  is  not  changed 
(off  axis  collecting  lens  moved  closer  to  transmitting  lens 
with  smaller  diameter  but  same  f/no.)  The  effects  are  essen- 
tially the  same  for  the  x^  distributions  as  above  in  B.4.2.3. 
(However,  the  fringe  number  remains  the  same  so  system  band- 
width increases) , 

B.4,2.5  Wavelength 

It  would  be  very  helpful  if  we  could  measure  r^(x)  at 
one  wavelength  and  determine  by  scaling  what  the  effects  of 
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switching  to  a laser  of  a different  wavelength  would  be. 

If  all  the  scatterers  of  interest  were  much  less  than  the 
optical  wavelength  in  diameter,  this  could  be  done  because 
the  scattering  intensity  from  all  Rayleigh  sized  scatterers 

is  proportional  to  x (the  product  of  i(a)X  ‘ is  propor- 

4 0 

tional  to  XQ  in  equation  B.2).  The  scaling  would  have  to 

include  a change  in  wQ  as  well  as  all  the  other  dependent 
factors  in  equation  B.2.  In  all  applications  of  interest 
that  seem  relevant,  however,  the  submicron  and  even  larger 
Mie  scatterers  contribute  significantly  and  the  distribu- 
tions are  not  scalable. 

B . 4 . 2 . 6 Optical  Geometry  Changes 

Due  to  the  complexity  of  the  probe  volume  geometry 
effects  and  the  Mie  scattering  functions,  most  changes  in 
the  optical  system  such  as  polarization,  scattering  angler, 
spatial  filter  size  or  shape,  etc.  do  not  produce  scalable 
changes  in  the  distribution  functions.  This  means  that  a 
thorough  parametric  study  of  the  behavior  of  the  functions 
we  have  been  describing  will  require  extensive  computer 
calculations  or  experimental  measurements.  Our  experiments 
have  demonstrated  that  the  power  law  behavior  of  the  dis- 
tributions of  X is  significantly  affected  by  the  optical  sys- 
tem parameters. 
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LIST  OP  SYMBOLS 


particle  radius 

PMT  anode  voltage  at  pedaBtal  peak^-pafticle 
radius  (irt  Appendix  B) 


accumulator  value  at • 

B bandwidth 

c »p<afcA«r''filter  aperture  width  - 2cwQ 

Cii(Vt2r*utocovariance  of  i(t) 

diameter  of  collecting  objective 
e electronic  charge 

f(  ) optical  probe  volume  response  function 

F focal  length  of  collecting  objective;  PMT 

collection  efficiency 

g^  random  single-photoelectron  charge  gain 

Planck's  constant 

impulse  response  of  PMT  and  succeeding  filters 
Fourier  transform  of  h(t) 
counting  index 

photocurrent  at  anode  or  succeeding  filter  output 
optical  intensity 
counting  index 

counting  index,  or  arbitrary  constant  in  scaling 
probe  volume  length 

time  varying  statistical  mean  current  <i  (t) > 

nknk-p  - nknk-q  or  other  function  of  n^ 
photon  correlation  sum  at  delay  pAx 
accumulation  at  (pAx,  qAx)  for  Dual  Correlate  Mode 
photoelectron  count  in  interval  (t^,t2) 
differential  particle  number  density 
photon  count  in  Ax  interval  about  kAx 
number  of  time  increments  Ax  (T  - NAx) 
number  of  particles  with  radius  a > a^  per  unit 


h 

h (t) 
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i (t) 
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j 

k 
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kp 

V 
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volumes 
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H 


I 


PX(X) 

px(x|y.2) 

P(A) 


probability  density  for  random  variable  x 

conditional  probability  density  for  x given  the 
values  of  the  random  variables  y and  z 

probability  of  the  event  A 

probability  that  the  random  variable  x exceeds 
the  value  x^ 


P(t) 

q 

r 

r, (t,  A) 
R<t) 

Rx(t,A>A') 

R,  (t) 

.it, 

t 

T 

T 

m 

u ( t) 

U(t) 

u 

V 

wo 

Wx 

X 

y 


z (t) 


optical  power  incident  on  photocathode 
integer  value  of  delay)  power  law  a^,  q real 
power  law  ar 

random  translational  location  of  jth  particle 
WRT  probe  center 

differential  amplitude  burst  rate 

instantaneous  rate  of  signal  bursts 

mean  rate  of  bursts  at  time  t with  A . > the  value 
00  J 

A'  given  by  / r^(A)dA  (in  Appendix  only)  t 

autocorrelation  of  A (t)  - <\(t)X(t+T)> 
a generalized  Poisson  shot  noise  signal  (Appendix  A) 
time 

total  data  collection  time 
inverse  of  mean  burst  frequency 
time  varying  velocity  component 
total  velocity  component 
mean  velocity  component 
vector  velocity 

l/e  intensity  radius  at  beam  waist 
effective  probe  volume  width  in  direction  of  flow 
position  variable 
positicn  variable 

random  y displacement  of  jth  particle  from  probe 
center 

random  probe  volume  amplitude  multiplier  due  to  y^ 
position  variable 
Poisson  impulse  process 
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GREEK 


a 


0 

«<t) 

Yj 

e 

Mt) 


n 


aj(t) 

T 

Th 

Tj 

e 


d> 

y 

w 


Uj 

id. 

m 

ck(X) 


integration  dummy  variable;  also  1/e  half-width 
of  burst 

integration  dummy  variable 

Dirac  delta  or  unit  impulse  function 

peak  value  X,  which  would  result  if  particle 
went  through-’the  probe  volume  center 

rms  error  ,**•*'' 

time  varying  mean  rate  of  photoelectron  pulses 

X'Tt)  - nP(t)/hv 

steady  background  photoelectron  rate 

peak  photoelectron  rate  (pedestal)  from  the  jth 
scatterer.  X^  » Yj*j 

optical  wavelength 

mean  signal  photoelectron  rate  ^ <X.>(see  equation 
3.24)  3 

product  of  photocathode  quantum  efficiency  and 
dynode  collection  efficiency 

optical  frequency 

2 

time  varying  variance  of  i(t);  <(i(t)-<i(t)>)  > 
delay  variable 

rise  time  or  pulse  width  of  low  pass  h(t) 

occurrence  time  for  jth  signal  burst 

beam  intersection  angle;  spherical  declination 
angle  from  direction  of  incident  light 

polar  angle 

Poisson  parameter 

Fourier  transform  variable  or  radian  frequency 
random  frequency  of  jth  burst 
mean  of  random  variable  uij 

normalized  differential  rate  function  r^(X)/|v| 


Special  Notation 

<x>  statistical  expectation  of  x 

00 

f(t)*g(t)  convolution;  / f (a) g (t-a) da 


00 

/ xpx  (x)dx 
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